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Abstract 

We present a description of calculations of the amplitude for e + e~ — > 
// process with account of electroweak and QED one-loop corrections. 
This study is performed within the framework of the project SANC. The 
calculations are done within the OMS (on-mass-shell) renormalization 
scheme in two gauges: in R^, which allows an explicit control of gauge 
invariance by examining cancellation of gauge parameters and search 
for gauge-invariant subsets of diagrams, and in the unitary gauge as a 
cross-check. The formulae we derived are realized in two independent 
FORTRAN codes, eef f Lib, which was written in an old fashioned way, i.e. 
manually, and another one, created automatically with an aid of s2n.f 
(symbols to numbers) software — a part of SANC system. We present a 
comprehensive comparison with the results of the well-known program 
ZFITTER for all the light fermion production channels, as well as with 
the results existing in the world literature for the process e + e~ — » ti. 
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Introduction 



The process e + e~ — > // with taking into account of the higher-order Standard 
Model (SM) corrections has been studied already several decades, yet in pre- 
LEP times. We mention papers from late seventies, PQ, 0, ED and @], which 
eventually lead to producing of precision theoretical predictions realized in 
the well-known computer codes of LEP-era: ZFITTER 0, BHM ALIBABA [7] 
and TOPAZO 8 . These codes played a very important role throughout the full 
life-time of LEP providing experimental community with precision tools for 
fitting of experimental data to the SM predictions, see, for instance 
A comprehensive review of underlying theory and methods which has been 
used for creation of these codes may be found in the monograph 11 . 

At LEP, the higher-order corrections for process e + e _ — > // might be 
studied ignoring external fermion masses, since the cms energy was far below 
it threshold. However, studies of finite mass effects has been started more 
than ten years ago |12j . [13] and |14| in connection with experiments at future 
linear colliders (see, for instance, the reviews [13] and |16|1- 

New wave of papers on finite mass effect was triggered by studies of the 
electroweak radiative corrections (EWRC) in the MSSM which may involve 
additional heavy fermions and sfermions, |17j-jl8j. 

The process e + e~ — ► it in reality is a six-fermion process (see JH]); one of 
the channels is shown in Fig. ^ 




Figure 1: The six-fermion e + e — * tt process. 

However, the cross-section of a hard subprocess, a(e + e~ — > it), with tops 
on the mass shell is an ingredient within various approaches, such as DPA 20 
or the so-called Modified Perturbation Theory (MPT), see [2"T] . 

One should emphasise that the treatment of even one-loop corrections with 
taking into account of finite-state fermion masses is extremely cumbersome 
and practically undouble 'manually'. Nowadays, all the calculations of such a 
kind are being performed with the aid of automatic computer systems, among 
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which one should mention first of all FeynArts package, which exists already 
more than ten years [22] and which is able to compute the cross-section of the 
process cr(e + e~ — > tt) at the one-loop level. 

In this article, we present a review of a new calculation (based on preprints 
E3 _ E3) 01 the e + e~ — ► // process at the one-loop level made with an aid of 
the computer system SANC where all the calculations from the Lagrangians up 
to numbers are going to be eventually automatized. This system is being cre- 
ated at the site brg . j inr . ru. It roots back to dozens of supporting form 26 
codes written by authors of the book JT] while working on it. Later on, the 
idea came up to collect, order, unify and upgrade these codes up to the level 
of a 'computer system'. Its first phase is described in [2"S] . 

One of main goals of this new calculation was to create a platform for 
the treatment of any 2/ — > 2f process within the SANC and to cross-check its 
results against the results obtained with the other existing codes for a rather 
complicated 2/ — > 2/ process at one-loop level taking into account final state 
masses. 

At present, SANC, like FeynArts, uses the OMS renormalization scheme, 
a complete presentation of which was made in jTJ. However, for the first 
time the calculations are performed in two gauges simultaneously: i?j and 
the unitary gauge. 

Note, that there was wide experience of calculations in the gauge for 
processes such as H — ► //, WW, ZZ, 7Z, 77, or e + e~ — * ZH, WW. So, in 
[23 and [2E] a complete set of one-loop counterterms for the SM is given. 
Electromagnetic form factors for arbitrary £ are discussed in [21] and . 
Explicit expressions can be found in the CERN library program EEWW [3T] . 

However, we are not aware of the existence of calculations in the gauge 
for the e + e~ — > tt process, although there are many studies in the £ = 1 gauge 

Q3-Q1 and C3-[IH] 

Additional purposes of this review are: 

— to explicitly demonstrate gauge invariance in R^ by examining cancel- 
lation of gauge parameters for gauge-invariant subsets of diagrams; 

— to offer a possibility to compare the results with those in the unitary 
gauge, as a cross-check; 

— to present a self-contained list of results for one-loop amplitude in terms 
of Passarino-Veltman functions Aq, Bq, Cq and Dq (as well as special 
functions eto, 60 1 Co and do which originate because of a particular form 
of the photonic propagator in the i?j gauge) and their combinations 
in the spirit of the book where the process e + e~ — > tt was not 
covered; this review may thus be considered as an Annex to the book, 
which completes its pedagogical aims for 2/ — > 2f processes at one loop; 
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— to provide a FORTRAN code, eef f Lib [32] for the calculation of the cross- 
section of this process for a complementary use within the MPT frame- 
work; 

— to compare the results derived with the aid of eef f Lib with the results of 
another code, which was created automatically using the s2n. f software 

— a part of SANC system, thereby benchmarking s2n . f software. 

This review consists of seven sections. 

In Section 1, we present the Born amplitude of the process, basically to in- 
troduce our notation and then define the basis in which the one-loop amplitude 
was calculated. We explain the splitting between QED and EW corrections 
and between 'dressed' 7 and Z exchanges. 

Section 2 contains explicit expressions for all the building blocks: self- 
energies, vertices and boxes, both QED and EW. Note that no diagram was 
computed by hand. They all were supplied by the SANC system. 

In Section 3, we describe the procedure of construction of the scalar form 
factors of the one-loop amplitudes out of the building blocks. One of the 
aims of this section was to create a frame for a subsequent realization of the 
renormalization procedure within the SANC project. 

Section 4 contains explicit expressions for the improved Born approxima- 
tion (IBA) cross-section and the explicit expressions for helicity amplitudes 
made of the scalar form factors at the one-loop level. 

Section 5 contains some additional expressions for different QED contri- 
butions that might be derived analytically. They are not in the main stream 
of our approach: Lagrangian — * scalar form factors — > helicity amplitudes — > 
one-loop differential cross-section. However, they are useful for the pedagogi- 
cal reasons, and their coding in complementary FORTRAN branches of eef f Lib 
provided us with the powerful internal cross-checks of our codes for numerical 
calculations. In reality, the eef f Lib version of February 2002 has three QED 
branches. 

Finally, in Section 6 we present results of a comprehensive numerical com- 
parison between eef f Lib and ZFITTER. Here we also present a comparison 
with another of our codes, which was created automatically using the s2n.f 
software. We also present a comprehensive comparison between the results 
derived with our two codes and the results existing in the world literature. In 
particular, we found a high-precision agreement with FeynArts results up to 
11 digits for the differential cross-sections with virtual corrections, and with 
recent results of j^H] within 8 digits, even when the soft photon radiation is 
included, see p^fj-pTH]. 
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1 Amplitudes 



1.1 Born amplitudes 

We begin with the Born amplitudes for the process e+(p + )e~(p_) — > 
f(q~ )/(<?+), which is described by the two Feynman diagrams with 7 and 
Z exchange. The Born amplitudes are: 



-4 

At = 



(1.1) 



7m 
1 



'7k 



/[ 3) 7+ - 2Q f s 2 w 



Ml 



+ I i 3)S fl^l+ ® 7/i + <W/7a« ® 7^ 



4 3) 4 3) 7m7+ ® 7^7+ + <UfV ® 7^7+ 



(1.2) 



where 7 ± = 1 ± 75 and the symbol ® is used in the following short-hand 
notation: 

7 M (£i7+ + Qi) O 7 „ (^27+ + Q2) 

= v (P+) 7m (^17+ + Qi) u [p-)u (q-)j v (i 2 7+ + Q2) v (q+) ; (1.3) 



furthermore 



(1.4) 



Introducing the LL, QL, LQ, and QQ structures, correspondingly (see last 
Eq. p. 211 ). we have four structures to which the complete Born amplitude may 
be reduced: one for the 7 exchange amplitude and four for the Z exchange 
amplitude ( 7/J <g) 7/J structure presents in both amplitudes). 



1.2 One-loop amplitude for e + e — > ff 

For the e + e~ — > // process at one loop, it is possible to consider a gauge- 
invariant subset of electromagnetic corrections separately: QED vertices, AA 
and ZA boxes. Together with QED bremsstrahlung diagrams, it is free of 
infrared divergences. The total electroweak amplitude is a sum of 'dressed' 7 
and Z exchange amplitudes, plus the contribution from the weak box diagrams 
(WW and ZZ boxes). 

Contrary to the Born amplitude, the one-loop amplitude may be paramet- 
rized by 6 form factors, a number equal to the number of independent helicity 
amplitudes for this process if the electron mass is ignored and the unpolarized 
case is studied 1 . 

1 If the initial-state masses were not ignored too, we would have ten independent helicity 
amplitudes, ten structures and ten scalar form factors. 
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We work in the so-called LQD basis, in which the amplitude may be 
schematically represented as: 

[i^+Ft (a) + i ltl F e Q (s)] ® [iw + F[ (a) + i^F* (a) + mjID^f (a)] ,(1.5) 
with 



D„ = (g+ -?-)„. (1.6) 
Every form factor in the Re gauge could be represented as a sum of two terms: 

FUo (s) = Fgl B (a) + F^ D (a) . (1.7) 

The first term corresponds to the £ = 1 gauge and the second contains all £ 
dependences and vanishes for £ = 1 by construction. 

The LQD basis was found to be particularly convenient to explicitly 
demonstrate the cancellation of all ^-dependent terms. We checked the can- 
cellation of these terms in several groups of diagrams separately: the so-called 
A, Z, and H clusters, defined below; the W cluster together with the self- 
energies and the WW box; and the AA, ZA and ZZ boxes. Therefore, for 
our process we found seven separately gauge- invariant subgroups of diagrams: 
three in the QED sector, and four in the EW sector. 

The 'dressed' 7 exchange amplitude is 

A™ = i^QLa(a)^® lfi , (1.8) 

which is identical to the Born amplitude of Eq. i|1.2fl modulo the replacement 
of a(0) by the running electromagnetic coupling a(s): 

OL 

a ( s ) = en-, ; T • (1-9) 



a 

1 

47T 



n£00-n£(o) 



In the LQD basis the Z exchange amplitude has the following Born-like 
structure in terms of six (LL, QL, LQ, QQ, LD and QD) form factors: 

AT =1 e 2 ^fl | 7m7+ ® w + F LL (a, t) + 7m ® w+F QL (a, t) 

+ 7^7+ ® T/j F L q (a, t) + 7 „ ® ltl F QQ {a, t) (1.10) 

+ 7^7+ ® (-imfDf,) F LD (s, t) + 7m ® {-imfD^) F QD (s, t) I , 
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where we introduce the notation for F- (s,t) : 



F,, 


(s,t) 




Fql 


(s,t) 


= 5J { f 3) F QL (s,t), 




(s,t) 


= I^6 f F LQ (s,t), 


F 

QQ 


(*»*) 


= S e S f F QQ {S, t) , 




(*»*) 


= iPlf ] F LD ( S ,t), 


p 

± QD 


(s,t) 


= S e lf } F QD ( s ,t). 



(1.11) 

Note that tilded form factors absorb couplings, which leads to a compactifi- 
cation of formulae for the amplitude and IBA cross-section, while explicit ex- 
pressions will be given for untilded quantities. The representation of Eq. (|1.10|) 
is very convenient for the subsequent discussion of one-loop amplitudes. 

Furthermore, in Eq. (|1.1U|) we use the Z/7 propagator ratio with an s- 
dependent (or constant) Z width: 

Xz{s) = if* — sr tt~ ■ (L12) 



2 Building Blocks in the OMS Approach 

We start our discussion by presenting various building blocks, used to con- 
struct the one-loop form factors of the processes e + e~ — ► // in terms of 
the Aq, Bq, Cq and Dq functions. They are shown in order of increasing 
complexity: self-energies, vertices, and boxes. 



2.1 Bosonic self-energies 

2.1.1 Z, 7 bosonic self-energies and Z— 7 transition 

In the i?£ gauge there are 14 diagrams that contribute to the total Z and 7 
bosonic self-energies and to the Z-7 transition. They are shown in Fig. [2J 

With S zz , S Z1 and S 71 standing for the sum of all diagrams, depicted by 
a grey circle in Fig. El we define the three corresponding self-energy functions 

S zz = (2tt) 4 z— 4^E ZZ , (2.1) 
S Z1 = (27r) 4 z-f^-£ Z7 , (2.2) 

107T C\y 

4 9 2s w 

£77 = (2tt) jS 77 . (2-3) 
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/ 

Z{l) ^ Z{j) 

VWW WVTO = E/| 

it v_y v 




(1) / 




+ VWV* WWW + 



WWW www 



(2) (3) ff 



+ vwv# www + 



WW*, www 

(4) 0- (5) W~ 



+ ww# www + ww* www 

(6) H (7) 0" 

.•*"•. 

+ VWV* #AVW + WW* <jMA/W 
*•--*'' '•-**' 

(8) X- (9) X+ 

(10) W (11) if 




+ WWWWWWVW + WWWWWVWW 
(12) 0+ (13) 00 

/ \ f \ 

II II 

+ vwvwW/www + wwwwwwwv 



( 14 ) A (2) 

+ wwwwwwvw 

Figure 2: ^(7) boson self-energy; Z-7 transition. 
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All bosonic self-energies and transitions may be naturally split into boso- 
nic and fermionic components. 

• Bosonic components of Z, 7 self-energies and Z-7 transitions (see dia- 
grams Fig. [2J) 



3 00 



9c! 



(l + 2ct)^--i- 



3 Cv 



1 

— f* r 
2 



i + S 



(2.4) 

bos ' F (s), 



->bos,F 



M 2 

(s) = ^ 
v ; 12 



44 (5 - 8c 2 „ - 124) + (1 - 44 ~ 364) 



xBS(-s;M w ,M, 
1 

' 18 



10 - 2r HZ + {r HZ - if R z 



Bg(-s;M H ,M z ) 



+2c 2 



+ 1 + (1 - r H1 

Z 

7-{l-r HZ )R 
18 



+ 1 + 84 - 244 L^Mi) 



4 (1-24) ^--6 (1 + 24)-^- 

-3(l + 24)-9r ffz -(l-r ifz ) 2 i? 2 



Here L^M 2 ) denotes the log containing the 't Hooft scale [i: 



L^M 2 ) = In _ ., 



M 2 
V 2 



(2.5) 



(2.6) 



and it should be understood that, contrary to the one used in |llj . we define 
here 



B ( - s; Mi, M 2 ) = J + < ( - s; M t , M 2 ) , 



(2.7) 



meaning that £? F also depends on the scale /j,. We will not explicitly main- 
tain n in the arguments list of and Bq . Leaving fj, unfixed, we retain an 
opportunity to control /i independence (and therefore UV finitencss) in nu- 
merical realization of one-loop form factors, providing thereby an additional 
cross-check. 
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Next, it is convenient to introduce the dimensionless quantities IT^° s (s) 
and lTr!° s (s) (vacuum polarizations): 



s|° B ( s ) = - s n^( s ), 



(2.8) 

£*° S ( S ) = sU^ s ( S ). (2.9) 
In Eqs. 1)2. 5fl and (|2.7|l and below, the following abbreviations are used: 



M 2 m 2 M 2 M 2 

lvl w _ "H p _ lvl w c, _ lvl z (0 1n \ 

Mi mi s s 



Since only finite parts will contribute to resulting expressions for physical 
amplitudes, which should be free from ultraviolet poles, it is convenient to 
split every divergent function into singular and finite parts: 



lfi° s ( 



3 _ + n hos < F 



- (*)> (2.11) 
t 

(3 + AR XV ) B F ( - s; M W ,M W ) + 4R w L^Mi), (2.12) 



and 



lC s (s) = I ~ + 34 + 2R W I ^ + n 



n 



bos,_F 
z 7 



(*) = 



1 



g+3r- +-1 | -+r 
1 

"9 



-"W ] ^ w 



\ - Ac 2 v ) R W L^(M, 



*° s > F (s), (2.13) 
B F (-s;M w ,M w ) 
). (2.14) 



With the Z boson self-energy, S zz , we construct a useful ratio: 

1 S zz (,s) - E zz (Af|) 



2>z (*) 



M| - a 



(2.15) 



which also has bosonic and fermionic parts. The bosonic component is: 

1 



V h z ° s (s) 
with the finite part: 



1 1 

6 + 3 



cl + 3ci 



V 



hoa.F 



(s), (2.16) 



V 



bos : F 



00 



1 4 2 17 4 ,6 

— c c — Ar 

12 3 w 3 w w 



(2.17) 



Ml 
M 2 - 



B F ( - s; M W ,M W ) -B F (- M 2 Z - M W ,M W ) 



12 



1 1 



777 - ~4 - 3ct <( ~ s; M w , M w ) 



12 

1 r HZ H r z 

3 22 «z 



Ml 
1 

12 



( - s; M H , M z ) — Bq( — M 2 Z ; M H , M z 



1 - (1 - r HZ 



Bg(-s;M H ,M z 



-R z (1 - r HZ ) [r HZ (L„(Af|) - l) - L„(Mf ) + 1 



• Fermionic components of the Z and 7 bosonic self-energies and of the 
Z-j transition 

These are represented as sums over all fermions of the theory, Ylf They, of 
course, depend on vector and axial couplings of fermions to the Z boson, vt 
and cif, and to the photon, electric charge eQf, as well as on the colour factor 
Cf and fermion mass mj. The couplings are defined as usual: 



v f 



(3) 

with weak isospin 1 j , and 



r(3) o^i.„2 



a f 



r(3) 



(2.18) 



Qf = —1 for leptons, 
Cf — 1 for leptons, 3 for quarks. 

The three main self-energy functions are: 



— for up quarks, — — for down quarks, 



(2.19) 



(y 2 - + a 2 f) sBf ( — s;rrif, m f^)+2ajmj B ( — s; mj , my) 



(2.20) 
(2.21) 

(2.22) 



The quantities rfe and II^ are different according to different couplings, 
but proportional to one function Bf (see Eq. (5.252) of m for its definition): 



H^( S ) = 4Y,c f Q 2 f B f (- S ;m fl m f ), 
f 

n z 6 7( s ) = 2 J2 c f®f v f B f(~ s;m f ,m f ) 



(2.23) 
(2.24) 
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As usual, we subdivided them into singular and finite parts: 



(2.25) 



IE. 



-.fcr.F 



(*). 



In Eq. Q2.25J1 . A/ = 24 is the total number of fermions in the SM. We do not 
show explicit expressions for finite parts, marked with superscript F, because 
these might be trivially derived from Eq. I|2.20|l and Eqs. I|2.23l) - I|2.24[) by 
replacing complete expressions for Bf and Bq with their finite parts Bj and 
Bq , correspondingly. 



2.1.2 W boson self-energy 

Next we consider the W boson self-energy, which is described, in the R^ gauge, 
by 17 diagrams, shown in Fig. 

First, we present an explicit expression for its bosonic component: 



•(*) = Ml 



^bos.F 



19 1 1 
~6RZ ~ 4 

(*)■ 



G 



2 I - — + 10 + 3r HW 



1 

e 

(2.26) 



where 



^bos.F 



(s) = 



Ml 
12 



l-40cl 



R w 



2 - 27 - Set, 



^ ^ + 8 i? w 



Bg(-s;M w ,M z 



-L + 2(5-r HW ) + (l~r Hn ,) 2 R v 



B^(-s;M w ,M H ) 



-8ai 



1 



— h 2 - i? w ) B 



;M W ,0) 



7-(l-r w )J2 w L M (Mi) 



18 1 
18 



1 - 164 + 4 I -5- +8 I #w 



- 7 I - [ — - 2 + r 



HW I 



4 1 
3 R w 
1 



1 



31 A +2 



9r\ 



3- + 6s w 3- - r HW (2 - r ffw ) 



(2.27) 
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u W+ 

(1) d (2) Z 

W+ W+ 

^1 ^% 

+ VWV| MAMA/ + VWV» %AM/V 

(3) 7 (4) ff 

0+ 0+ 
, »- ^ , »- ^ 

+ WW* #A/WV + VWvW AA/WV 

(5) Z (6) 7 

0+ 0+ 
, *- , 

+ ww» Wvwv + vwvW Wvwv 



\ / 



(7) (8) 

Y z,l X + 

.•♦*•-. 

+ vwv* *a/wv + vw\a* Wvwv 




+ VW\AAA^WWWV 
(16) ^ 



I I 



+ VWVW^WWWV + VW\AAA^/WWWV 
( 17 ) A 

+ vwvwwwwwv 

Figure 3: M 7 " boson self-energy. 
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Secondly, we give its fermionic component: 

s ww («) = - s E °f B f ( - s ; m f > m /) + E c / m / Bi ( - s ; m /' > m /) > ( 2 - 28 ) 

/=d / 
where summation in the first term extends to all doublets of the SM. 



2.1.3 Bosonic self-energies and counterterms 

Bosonic self-energies and transitions enter one-loop amplitudes either directly 
through the functions T> z (s), IL„(s) and n z _,(s), or by means of bosonic 
counterterms, which arc made of self-energy functions at zero argument, owing 
to electric charge renormalization, or at p 2 — —M 2 , that is on a mass shell, 
owing to on-mass-shell renormalization (OMS scheme). 

• Electric charge renormalization 

The electric charge renormalization introduces the quantity z 7 — 1: 

2 



n 77 (o) 



Mi 



S 3 q(0) 



with bosonic (see Eq. (6.161) of [lip: 



bos 



and fermionic 



(z 7 - 1) 



fcr 



|E^)i + n ^(0) 



(2.29) 



(2.30) 



(2.31) 



components. 

• p-parameter 

Finally, two self-energy functions enter Veltman's parameter Ap, a gauge- 
invariant combination of self-energies, which naturally appears in the one-loop 
calculations: 



v (m£)-e zz (m|) 



with individual components where we explicitly show the pole parts: 



bos 



Ap 
Ap iei 



41 



1 

6c^ 6 
3 cl, 2 



7c z 



1 



Ap 



bos.F 



Nf-isi^CfQ^l + Ap^. 



(2.32) 

(2.33) 
(2.34) 
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The finite part of Ap bos is given explicitly by 
1 4 17 



Ap 



bos.F 



12c$ 



4c 2 



£? F ( - Ml; M W ,M Z ) - c 2 w B$( - M 2 Z ; M W ,M W ) 
+ (l - \r HW + ^r 2 HW ) B F ( - ; M w , M H ) 

-Bg(-M 2 Z ;M Z ,M H ) 



4s 2 X(-JV4;M ff ,0) 



1 

12 



-24 + r w )i P (M 2 ) 



-3- + 14 + 164 - 48ct + r„ w ) L M (M^) 



i 2 2 

+ S w r HW 



L»{M 2 H ) - 1 



19 22 

3^ + y 



(2.35) 



while the finite part of Ap fcr is not shown, since it is trivially derived from the 
defining equation 1)2. 32[) by replacing the total self-energies with their finite 
parts. 

2.2 Fermionic self-energies 
2.2.1 Fermionic self-energy diagrams 

The total self-energy function of a fermion in the gauge is described by 
the six diagrams of Fig. 




/ /*\ f 




(1) 7 



(2) 



(3) W 



r 



(4) H 



(5) 0° 



(6) 



Figure 4: Fermionic self-energy diagrams. 
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A detailed calculation of the wave-function renormalization factors ^jz L R 
associated with these diagrams may be found, for instance, in [5]. 
It is convenient to distinguish the electromagnetic components, 



2s 



and the weak components, 

I = (w v ± W a ) 



where five non-zero contributions are: 



1 1 



In- 



(2.36) 



(2.37) 



1 1 



(vj + a} + 2a)r fz ) - + (vj + a}) — Bq ( — m 2 f ;M z ,m f ) 



r fz 



+ L M (M|)-1 + 2(l + 2r fz )M 2 z B 0p (-m 2 -m f ,M z ) - L^m 2 ) 



+ 2a 2 f r fz [—\B^(~m};M z ,m f )+L^M 2 z ) - 1 



6M 2 B { 
1 f 1 



■ 0p (-m 2 f ;m f ,M z ) - L^m}) + l\ i, 



(2.38) 



w, 



16 e 



2 + r 



(2 + r+) L^M 2 v )-r f/w L^m 2 ,) 



+ (2 + 2r fw +r- +r+r-) B^ ( - mj; M w ,m 2 fr ) -1 



+ 2M 



2-r+-(r-) 2 B 0p (-m 2 f ;M w ,m 2 f ,)\, 



ff/ [Bg ( - m 2 f] M H ,m f ) + ^ (M 2 } _ X ] 

2(4r /ff -l)MlB 0p (-rn^m / ,M H ) -L M (m^) + l|, 



4c?, 



1 1 



( - mj; M z ,m/) + L M (A/ Z 2 ) - l] 



+ 2B^( - m); M z , m f ) + L^m}) -2J, 



(2.39) 



(2.40) 



(2.41) 



Tfw 



(2-r") L^MjJ,) - rf W Lfj,(m 2 ,) 



+ (2 - 2r /w - r+ + r+r") ( - m 2 f ;M w ,m 2 f ) - 1 



(2.42) 



with = rt w ± r/' W and /' being the weak isospin partner of fermion /. 
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2.3 The Zff and 7// vertices 



Consider now the sum of all vertices and corresponding counterterms whose 
contribution originates from the fermionic self-energy diagrams of Fig. 01 This 
sum is shown in Fig. [5] 






Figure 5: Zff and 7// vertices with fermionic counterterms. 
The formulae which determine the counterterms are: 
i^ ct = 2(^-1), 



F2' ct 

F'Z.Ct 
c 

F, 



Q 

z.ct 



{^:-i)-{^r R -i), 

aj(^-I)-Sj(^-I), 



where 



v f ~ a /> 



07 = Vf+a f . 



(2.43) 
(2.44) 
(2.45) 
(2.46) 

(2.47) 



For the sum of all 7 — > // and Z — > // vertices (the toiaZ j(Z)ff vertex 
depicted by a grey circle in Fig. |5J) we use the standard normalization 



16tt 2 ' 



and define 



while we denote the individual vertices as follows: 

G^s) = G*(s) + G*(s) + G™(s) + G»(s), 
Z^{s) = Z*{s) + Z*{s) + Z™{s)+Z«{ S ). 



(2.48) 

(2.49) 
(2.50) 



(2.51) 
(2.52) 



All vertices have three components in our LQD basis. 
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2.3.1 Scalar form factors 

Now we construct the 24 = (4 : B = A, Z, H, W -virtual) ® (3 : L, Q, D) ® (2 : 
7, Z -incoming) scalar form factors, originating from the diagrams of Fig. [S] 
They are derived from the following six equations — three projections for 7// 
vertices: 

F2 B (s) = -^- m {G;(s)[ig 3 lfll+ } + s w Q f F2' ct ), (2.53) 

Sw * f 

*T00 = ^{g^^ + ^Q/F^ 4 }, (2.54) 
F2 B (s) = -AgyG* (a) \g 3 m t ID^] , (2.55) 

and three projections for Zff vertices: 

ffW = %{^Wfe 3 7 M 7 + ] + — F Z L A, (2.56) 
I) I Cw ) 

F z Q B (s) = Tf{^( s )[^] + ^ Ct } ; (2-57) 
n B (s) = ^ZZWtfmtlDp]. (2.58) 

Here we introduce the symbol [. . . ] for the definition of the procedure of the 
projection of (s) and Z^ (s) to our basis. It has the same meaning as 
in form language namely, e.g. G^ (s) [ig 3 7 M 7+] means that only the 
coefficient of [«ff 3 7 M 7+] of the whole expression G^ (s) is taken (projected). 

The factors l/(QfS w ), 2/(s K ,/j 3 ' ) ) and 2c^/(s lv /l 3 ^) for 7// vertices, and 

the factors 2c w /Ij , 2c w /Sf and 2c w /Ij for Z// are due to the form factor 
definitions of Eq. I|l 

The total 7// and Zff form factors are sums over four bosonic contribu- 
tions B — A, Z, W, H . All 24 components of the total scalar form factors in 
the LQD basis look like: 

00 = F2$* (s) + Fl% )z D (s) + F2$£ (s) + F2$S (s) . (2.59) 

The quantities F2^o (s) originate from groups of diagrams, which we will 
call clusters. 

2.4 Library of form factors for Bff clusters 

Here we present a complete collection of scalar form factors F2^d ( s ) origi- 
nating from a vertex diagram with a virtual vector boson, contribution of a 
scalar partner of this vector boson, and relevant counterterms. 
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Actually three gauge-invariant subsets of diagrams of this kind, A, Z and 
H, appear in our calculation. They may be termed clusters, since they are 
natural building blocks of the complete scalar form factors, which are the aim 
of our calculation. Again, in the spirit of our presentation, we write down 
their pole and finite parts. 

The remaining vertices with virtual W and <p + , <f>~ with relevant countert- 
erms we also define as the W cluster. However, the latter diagrams do not 
form a gauge- invariant subset. 

We note, that F2 A (s) and F2" (s) are equal to zero. 

2.4.1 Library of QED form factors for Att clusters 

Up to the one-loop level, there are two diagrams that contribute to the A 
cluster, see Fig. El 

Separating out pole contributions 1 /e, we define finite (calligraphic) quan- 
tities. We note that, if a form factor F%* (s) has a pole, then the corresponding 
finite part Ja{s) is /i-dependent. 




Figure 6: A cluster. The fermionic self-energy diagram in brackets gives rise 
to the counter term contribution depicted by the solid cross. 



Since the scalar form factors of A cluster become UV-finite after wave 
function renormalization, for all form factors, which are also separately gauge- 
invariant, we have: 

F2%}£=?l$t (2.60) 

Individual components are: 
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„ „ 8m? 



3r 



8TT7 2 /" ( 3) 



A 3r (5/ 



Bq'( - s;m/,m/) - ( - m 2 f ;m f ,0) 
Bq( - s;m/,m/) - B(f ( - m 2 f ;m f ,0) 



D 

with 



if A 3r 



Bg(- s;m f ,m f ) - Bg( - m 2 f ;m f ,0) 



A 



3r 



4m f — s . 



(2.61) 



(2.62) 



2.4.2 Form factors of the Z cluster 

The diagrams shown in Fig. contribute to the Z cluster. 






Figure 7: Z cluster. The two fermionic self-energy diagrams in the second 
row give rise to the counterterm contribution depicted by the solid cross in 
the last diagram of the first row. 

Separating out pole contributions 1/e, we define finite quantities: 

F2 Z (s) =Tl z (a) , 



K 00 



1 1 

16|Q7^ 



FT (») =K z (*) • 



(2.63) 
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Here the finite parts are: 
TT (s) ~Q f v f ^2 ( 2 + ^) M 2 z C (-m 2 f ,-m 2 f ,-s;m f ,M z ,m f ) 

-3Bg(-8;m f ,m f ) + 2Bg (- m 2 f ;m f ,M z ) - L^mj) (2.64) 

M 2 1 
+ B$ 1 (-m 2 f ;m f ,M z )-2{l + 4r fz )-^-L ab (m f ,m f ,M z )y 

T V W = ^fi ( 4(1 ~ r i z) + ^) M ' C °( _ m 2 f ,-m 2 f ,-s;m f ,M z ,m f ) 
-3B^(-s;m f ,m f ) + AB^ (- m};m f , M z ) +L„(Af|)- | 
-^Sj 2 (-m^;m/,M z ) - 2 (1 + 2r /z ) M|B 0p ( - m^; m/, M z ) 



+ 2vfdfrf Z 



-4M 2 Ca( — m 2 , —m 2 , — s; rrif, M z ,m,f) 



+ 2 (L„(m?) - L„(M|)) - 2 (1 - r fz ) B£( - mj; m /; M z ) + 1 



»7* V 

2 



( (1 + 2r/ z ) Mf B 0p ( - m^; to/, M z 



+ 2a f rf z 



Bg(- s;m f ,m f ) + L„(M 2 ) - r fz B% 2 ( - m 2 f ;m t ,M z ) 



+ 6M 2 B Qp (-m 2 f ;m t ,M z ) 



M| 

-— iab (m/,m/,M z ) 
A 3r 



(2.65) 



1 J 7} 3) ^A 3r \ 2 



— 4MjCo( — To/, —to/, — s; m/, M z , m/) 



+ Bjf ( - s;m f ,m f ) - 2Bq ( - m 2 ;m f ,M z ) - L^(m 2 ) 



Ml 



+ B% 1 (-m 2 f ;m f ,M z ) +2 + 6— ^L ah (m/,m/,M z ) 



i3r 



+ 4 



2 3r /2 



M 2 C ( - m 2 , -m 2 , -s; to/, M z ,m f ) 



+ B^(-m 2 f ;m f ,M z )+L fl (M 2 z )-l 
~ r fz Bg(- s;m f ,m f ) + L^m))- 



2( 2-3^- )L a6 (m / ,m / ,M z ) 



(2.66) 
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4c4 



R z 



x MjCo( — mj, —mj, — s; to/, M z , m/) 
-9B^(-s;to/,to/) +8B^(-m 2 ,;TO/,M z ) -L fJi (m 2 f )-2 

+ B$ 1 (- m 2 ;m f ,M z ) - 2 (1 + 2r /z ) M| B 0p ( - to 2 ,; m/, M z ) 

2 r 

--a^ 4m 2 f Co(-m 2 f ,-m 2 f ,-s-,mf,M z ,mf) (2.67) 

+ 3r /z Bq (~ s;m f ,m f ) + L^rn 2 ) + Bq ( — to 2 ; to/, M z ) 
+ 3L^(M 2 ) + 2L^m 2 ) - 1 + 2B£( - to 2 -; to/, M z ) 
+ 2 (1 - 7r /z ) M 2 £ 0p ( - to^to/, M z ) 



- 2 



(3^ + a/) (1 + 4r /z ) - 2a 2 f r f2 



1 fl 



M 2 

-r^-L ab (m f ,m f ,M z ) \, 



2 (2-2r„ + -j M 2 C (-m 2 f ,-m 2 f ,-s;m f ,M z ,m f ) 

-3Bg(-8;m f ,m f ) +4Bg(-m};m f ,M z ) + L^{m)) -2 

- B£( - to 2 -; to/, M z ) - 2 (1 + 2r /z ) M 2 B 0p ( - m 2 f ;m fl M z ) 

M 2 

- 2—^-L ab (mf,mf,M z ) 

^3r 



+ «/ r /z [ Vf 
1 



r fz 



2M 2 Co( — rrif, —ntf, — s; to/, M z , to/) 
-B(f ( - to 2 ; m/ , M z ) + L^m 2 ) - 1 - - to 2 ,; to/, M z ) 



M 2 



(1 + 2r /z ) M 2 B 0p ( - to^to/, M z ) + 6— ^L a6 (to/, to/, M z ) 



i3r 



7 a / 



8M 2 Co( — to 2 , —to 2 , — s; m/, M z , to/) 



- s; to/, to/) - L^m}) + 2 + fl£ ( - to 2 ,; to/, M z ) 

M 2 

6M 2 B 0p ( -m 2 f ;m f ,M z ) - 10—^L ab (m f ,m f ,M z ) 



5/ 



4M 2 C ( 



-TO/,-s;to/,M z ,to/) 



M 2 1 \ 1 

-B^ , (-s;m / ,m / ) + 1 - 6^i o6 (m/, to/, M z ) , (2.68) 
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K z (s) 



2lf ] c 2 w A 3 r 



(34 + vj) Vf 



— 4M z Cq( — rnj, —nit, — s; to/, M z ,mf) 
+ Bg(-8imf,m f )-2Bg(-m};mf,M s )-L lt (m})+2 



Ml 



+ B£( - m};m f ,M z ) + &-^L ab (m f ,m f , M z ) 



3r 



+ 2v f a) 2 \7r fz ~ M 2 Z C ( - m 2 f , -m}, -s;m f , M z ,m f ) 

+ Bg(- m 2 f ;m f , M z ) + L^{M 2 Z ) - 1 



r fz 



Bq( - s;m f ,m f ) + L^mj) - 2 



2 4-3- 



i;ir 



Lab {mj,mf,M z 



(2.69) 



In Eq. I|2.65|l the 'once and twice subtracted' functions B^ and B^ 2 are met 



B, 



,n\ -n,y,m f ,M s ) =— \b${ - mj; m s ,M z ) + L^M* ) - ll , (2.70) 



1 



ro/,M,)+L M (Aff)-l 



r /z L /i (m^)-L M (M^) + - 



They remain finite in the limit to./ — > 

We note that, for the Z cl 
are separately gauge-invariant 



We note that, for the Z cluster, all the six scalar form factors F2,q,£ (s) 



2.4.3 Form factors of the H cluster 

The diagrams of Fig. |S] contribute to the H cluster. Separating UV poles, we 
have: 

*T (s) =TT (s), 



(2.71) 
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(z)h/\ {z)<f; 

t vwv»' \f + ww»' f/ 




-LTU- 

\ / 

H 

Figure 8: H cluster: the vertices and the counterterm. 

with the finite parts: 

Fq* (s) =^r fw ^8mfC ( - mj, -mj, -s; m f , M H ,m f ) 

+ B£(-s; m f ,m f ) + L^mj) - 2 - B£( - mj; m f ,M B ) 

-2(l-4r fH )MlB 0p (-m 2 f ;m f ,M H ) 

M 2 1 
-2-^L ab (m f ,m f ,M H )\, (2.72) 

H" (*) - - ^-{~6M 2 H Co( - mj, -mj, -s; m f , M H ,m f ) 

+ W£(- s;m/,m/) - 4B^( - m^; m/, M H ) - L^mj) 

M 2 1 
+ 2 + Bj 1 (-m^;m / ,M H ) +6^L a6 (m / ,m / ,M„)J, (2.73) 



26 



Fl" ( s ) =^ r fw^rn 2 C (- m 2 f , -m 2 f , -s;m/, M H ,mf) 



A(l-r fz ) + (l-r HZ ) 2 R z 

x M 2 C {-m 2 f ,-m 2 f ,-s-,M H ,m f ,M z ) 

+ 2B^(-s;M z ,M H ) - ^ { - s; m f , m f ) + ^L^m 2 ) + 2 
- \ B l\ ( - m 2 f ;m f , M H ) - (1 - Ar fH )M 2 H B 0p ( - m}; m f ,M H ) 

+ {l-r HZ )R z B^(-m 2 f ;M z ,m f ) - B^(-m 2 f ;m f ,M H ) 
Ml - 



+ 



Kir 



(r HZ - 8r/ z ) L ab (to/, m fl M H ) 



2 (3 - r HZ + Ar fz ) L Hi {m f , M H ,M Z ) 



}■ 



(2.74) 



+ 



■ff (s) =r /w |mjC ( - m 2 , -m 2 , -s;m f , M H ,m f ) 

1 + (1 - r HZ ) R z M 2 C ( - m), -m), -s; M H ,m f ,M z ) 

Bq{ ~ s;m f ,m f ) + L^(m 2 ) - 2 
Jb£(- m)- M z ,m f ) - Bg( - m 2 f ;m f , M H ) 



1 



+ R 

- - ; m/, M a ) - J (1 - 4r /H ) M 2 H B 0p { - m 2 ;m f , M H ) 
1 M 2 

- - A ir~ L o.b {m f ,m f ,M H ) 

+ \j- [y rfz + (3 + r HZ ) (1 - r HZ )i?. 

x M 2 Co(-m 2 ,-m 2 ,-s-,M H ,m f ,M z ) 
+ Bg(-s;m f ,m f ) -2B$(- s;M z ,M H ) -3 

+ (3 + r„ z ) i? z [b£ ( - m); M z ,m f )-Bg(- m 2 f ;m f ,M H ) 
„M 2 



i.'ir 



(1 - Arf H )r HZ L a b {mf, mf,M H ) 



- (3 - r HZ + Ar fz ) L m {m f ,M H ,M z ) 



)}• 



(2.75) 
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K H is) 



-3r fz M 2 H C { 



-m 2 ,-s; m/, M H ,mf) 



2(r H2 



-hz + 2r /z 

Vf 2 



x MjC (-m 2 f ,-m 2 f ,-s-,M H ,rn f ,M z ) 



Bg(-m 2 f ;M z ,m f ) - B^ (- m 2 f ;m f ,M H ) 

-2Bg(- s; M Z ,M H ) + 25 F ( - m^; M z , m/) 
3 r 

Bq( - s;m/,m/) - B F ( - m 2 f ;mf,M H ) 



+ 2 rHZ 



B* ( - m 2 ;m f , M H ) + L^(M 2 H ) - 1 
B$ ( - mj; m/, M H ) + L p (roJ) - 



M 2 

+ Srfz—^Lat (mf,mf,M Hi 

A 3r 



(2.76) 



Again, the five (one does not exist) scalar form factors in Eq. (|2.71|) are 
separately gauge- invariant. Note also that UV poles persisting in the scalar 
form factors of the H cluster cancel exactly the corresponding poles of the Z 
cluster. In other words, the form factors of the 'neutral sector' cluster (Z + H) 
are UV finite. 

In total, we have 11 separately gauge-invariant building blocks that origi- 
nate from Z and H clusters. 



2.4.4 Form factors of the W cluster 

Finally, the W cluster is made of the diagrams shown in Fig. |5J 

In the formulae below, we present contributions to scalar form factors 
from all the diagrams of the W cluster, not subdividing them into abelian 
and non- abelian contributions. To some extent two sub-clusters are automat- 
ically marked by the type of arguments of Co functions and typical coupling 
constants. Separating poles, we have: 



p-yw 


(«) 




Q 


(«) 




rp-yw 

D 


(*) 




K w 


(«) 






(«) 


=n w w > 




(*) 
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/' + 




+ 



/' + 




Figure 9: W cluster: the first row shows the abelian diagrams of the cluster, 
the last row the non-abelian diagrams; the second row shows diagrams that 
contribute to both counterterm crosses (last diagrams in first and third rows). 



Here the finite parts are: 
-i(6 + r 



■W J 



M^C (- m 2 ,-m 2 ,-s\ mf,M w ,m f ,) 



2 
1 



B o{~ s;mf>,m f >) - Bq( - m 2 f ;m f ',M w ) 
+ \{*-r-) [B^ w {-m 2 f ;m r ,M w )-l 



+ 



(2 + r ) (3 + 3r fw +r f>w ) 



M 2 1 
-r-^Lab {mf,mj>,M w ) 

^3r J 



m 2 ,—m 2 7 — s; M w ,mf> , M w ) 



-r f , w (2-r-)MlC ( - — 

+ \ (6 - r~) [fi F ( - s- M W ,M W ) -Bg(- m 2 f ;m f , , M w ) ^ 
+ 2B%{- m 2 ;m f , 7 M w ) + \ (2 - 3r") \b% w ( - m 2 f ;m r ,M w ) + 1 



- U- 



1 M 2 

2 + r fw (U + r fw )-r rw (l + r rw ) 

J A 3r / 



_Qs_ 

47 



x L na (m f ,m f ,,M w ) 

^(A^ Im - 3)zlm<( - m 2 ;M w ,m f ,), (2.78) 
'/ 
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where A^- Im = (1 - rf, w ) - — - 

The last term in Eq. I|2.78[l is due to a non-cancellation of the imaginary 
part of the function Bq ( — m 2 \ M w , my) which appear in real counterterms 
and complex- valued vertices. 

"Zq7 { Qf ' ~ 4M - C °( " m /' _m /' ~ s; ™/'' M «" m /0 

+ B,f ( - s;m/',m//) - ( - m 2 ;mf>,M w ) - 1 

_i±2V^ w (_ m 2; m/ ,,M w ) 

(l-r / ^)(2 + ^) _ 1 _ r/w + 2r//w 

x M 2 v B 0p ( - m 2 f ;m f ,,M w ) 

M w 

+ 2 (3 + r ) —^-L ab (mf,mf,M w ) 

&3r 



2/f 

2 + ?*/'»- „_f 



2m 2 ,Co( — my, — my, — s; M w , to/', M w ) 
-a;M w ,M w ,) + B^(-m 2 ;m f 
B dw( ~ m 2 f ;m f/ ,M w ) 



>,M W )-1 



Ml 



x M 2 v B 0p (-m 2 f ;m f ,,M w ) +2(7 + r+)^L„ Q (m/.m//, Af„r) 



3r 



-A^r Im zImB^( - m 2 ;M w ,m f ,) , 



(2.79) 



Qf 

A 3 , 1 2J} 3) 



8 + r~ 



6r// 



2 

Rw 



x Af^C ( - m^, -my, -s; m,f>,M Wi my/) 

+ (10 + r/ w - 3ry/ w -) Bq( - s;m, f ,,m f ,) - Bq ( - m 2 ;m,f>, M w ) 
+ (2 + r+) [Sj w ( - mj ; m fl , M w ) + 1 

+ 6 {l + r fw ){2 + r fw ) -r f , w (l + r f , w ) ^p^L ab (m f ,m f , , M w ) 
L J A 3r 

+ 2 (1 - r+) (2 - r fw ) - r f . w (l - 4ry, w - JL 



M 2 



30 



+ (6 - r fw - 5r f > w ) [b£ ( - s; M w , M w ) -B$(- m);m f ,,M w ) 
+ (2 + r+) \b% w (-m 2 f ;m r , M w ) - 1 



- 6 



(1 - r fw ) (2 + r fw ) - r f , w (1 + 2r fw + ry w ) 



M 2 

x -r-^ L ™a (m f ,m f >,M w ) 
A 3r 



(2.80) 



7' 



4/ 



(3) 



3+(l + r-) 2 + A 



M W C ( - m 2 f , -m 2 , -s; m,f,M w ,mf) 



--(6 + r ) Bg(- s;m f/ ,m f/ ) - (- m 2 f ;m f , , M w ) 



2-r" 



B lw ( - m /; m /' > M »-) - 1 



+ 



4- (2 + r~) (3 + 3r fw +r f/w ) 



M 2 

-r^Lab (mf,mp,M w ) 



+ r f , w MlC ( 



},-s;m f ,,M w ,m f/ ) 

+ ^r-Bg(-s;M w ,M w ) -r fw B* '(- m 2 f ;m r ,M w ) -1 

+ r// w Bo(-s;mf>,mf>)-2 - (2 + r f , w ) B% w (- m 2 ;rn f ,,M w ) 

- (! - r /w) (2 + r /w ) - r/- w (1 - 2r /vv + r/- w ) 

x MlB Qp ( - m 2 f ;m r , M w ) 

M 2 1 

- 2r/- w (l + r") -r-^£a& (m f ,m f/ ,M w ) } 

^3r J 



- c„ 



4(1 ~r fw ) + -r// w 



s 2 -r 2 
4 + w w 



(4 + r") 



x MlCo(-m 2 f ,-m 2 f ,-s;M w ,m f/ ,M w ) 



+ - (2 + r") < ( - m 2 ; M w , M w ) - B F ( - m 2 ;m f ,,M w ) 



--(2-r") |fl£ H ,(-roJ;m / /,M H ,)+l| - 2B F ( - to 2 ; m r , M w ) 
-\U+ 12r fw - ^-j^r fw (7 + r fw ) 



2 V 

rfw 



s 2 -c 2 
4+ £w^ (1 _ r// 



*3r 



■L no (mf,mf,M w ) 
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+ ^y^y (^<T f - 5 f Ajr iin ^jiImBv( - m 2 f ; M w ,mp) , 

- ^^^MlCo{-m},-m 2 f ,-s-,m r ,M w ,m r ) 

- Bq ( - s;mp,mf/) + Bq ( - m 2 ;m f >,M w ) + 1 

M 2 1 
-2(3 + r-)-^L o6 (m / ,m / ,,M w ) ^ 

^3r J 



(2.81) 



(2 + rp w ) B% w ( - m 2 ;m f < , M w ) 



(l-r fw )(2 + r-) +r f , w r~ 



MlB 0p (-m 2 f ;m f ,,M w ) 



+ ifcl^f ^ 2 ^ (r f >wMlC (- to 2 , -to 2 , -a; M w ,m f ,,M w 
B^(- S ;M w ,M w )-B^(-m 2 ;m r ,M w ) + l 



s 2 -c 2 



(3 + r fw + r f , 

w ) 



Ml 



+ jA?- Iln ilmBQ ( - m 2 f ;M w ,m fl ), 



^3r 



■L na (mf,m f >,M w ) 



(2.82) 



8+(r ) 2 -6r rw + ^- 
it w 



A 3 rl2/f 

x M^C ( - m 2 ,-m 2 , -s;m// , M w ,m { >) 
-^(W + r fw -3r f , w ) Bg(- s;m f ,,m f , 

-l(2 + r+) [ J B| w (-m2;m / ,,M w ) + l 



m 2 ;mp , M w ) 



(1 + (2 + r /vv ) - r// w (1 + r// w ) 
+ r/' W ^M^C ( - to 2 , - to 2 , -s; m f >,M w , to//) 



M 2 

— -L 6 (m f ,m fl ,M w 

^3r 



+ 



- s;m/',m/') - Bjf ( - m 2 f ;m f >,M w ) - 1 



M, 



- 2- B dw( - m 2 f ;m f ,,M w ) - 3 (l + r ) -^L ab (m f ,m f ,,M w ) 



- c„ 



4r// 



3 2 -C 2 



(1 - r/ w ) (2 - r/ w ) 
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■rf> w (5-r fw -4r f , w - jj- 



x M^C a (-mj,~m 2 f ,-s;M w ,m f ,,M w ) 



4- 



(4 - r fw - 5r//w) 



B^(-s; M w , M w ) -B$(- m};m f >,M w ) 
B dw(- m p m f'' M w) 



s 2 - c 2 

4 ,2 r tb 



Ml 



6(l-r+) -^L m {m f ,m f ,,M w ) 

^3r 



(2.83) 



Here we introduce more ratios, which were not given in Eq. (|2.1U|) : 

r ± =r fw ±r } . w . (2.84) 
Furthermore, we used one more 'subtracted' function: 

Bf w (-mj;m fl ,M w ) = — ((1 - r f , w ) \b% ( - m};m f ,,M w ) (2.85) 



+L„(M 2 ) - 1 



7'w 



L,(m 2 ,)~L,(Ml) 



and the three auxiliary functions: 



L«6(Mi, M 2 , Af 3 ) = (M| + M 2 - M 2 2 ) C„( - m 2 ,, -to 2 , - s; M 2 , M 3 , M 2 ) 

- < ( - s; M 2 , M 2 ) + < ( - to 2 ; M 2 , M 3 ) , (2.86) 
(M 3 2 - M 2 - M|) C ( - m 2 - , -to 2 , -s; M 3 , M 2 , M 3 ) 
+ < ( - s; M 3 , M 3 ) - < ( - to 2 ,; M 3 , M a ) , (2.87) 



L na (M l7 M 2 ,M 3 ) 
L m (Mi,M 2 ,M 3 ) 



- (M| + M 3 2 ) - 2M 2 



(2.88) 



x C ( - mj, -m 2 ,, -s; M 2 , Mi,M 3 ) + B$ ( - s; M 3 , M 2 ) 
- ifijf ( - m 2 ; Mi, M 2 ) - i< ( - m 2 ; M 3 , Mi). 

Four scalar form factors, !Fq% , as follows from calculations, are both gauge- 
invariant and finite, thus enlarging the number of gauge-invariant building 
blocks to 15. On the contrary, two form factors, T2^ w , are neither gauge- 
invariant nor finite. Gauge dependence on £, as well as UV poles, of L form 
factors cancel in the sum with the WW box and the self-energy contributions. 
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2.5 Library of scalar form factors in the limit rrif = 

In the limit mf — only nine scalar form factors of vertex clusters give non- 
zero contributions: Fq A , F^ A Q , F2 Z Q , F* z q and F2' zw , which, in turn, may 
be expressed in terms of three auxiliary functions F^ m , F]f m and F^ m . These 
nine form factors are: 



<31im 



± L lim 



Olim 



jpzz 
tlim 



P7W 
± L lim 



L lim 



p2A P ZA Cl z a Z P A 

— r L\im — r Q\i m ~ ^f b w r lim i 

-2Q f v f a f F£ m , 



- p 2Z _!, A 2 P z 

Qlim ~~ /i„2 "rl™ ' 
4t w 



1 



4/l 3) c 2 



(3vj + a}) a f F{( m . 



Qf F w 

(3) lim 



21 



r// w (2 + r/' W ) + 4d w + d^(2 + r f > w )R w 

f ' 



1 -Rw 



BS(-s;M w ,M w ) +L„(m 2 f ,)-l 

2i + i +1+(2+ ^ )i? " 



4/ 



2L„(m?,) 
Fl 



"Tf'w 



2d v 



(3) 



lim 



1 

~ 2 

-i? 



(2 + d 2 v R w )M 2 / C ( - mj, -m 2 f , -s;m,f,M w ,mf) 

B^(-s;M w ,M w )-B^(-s;m f ,m f )+2 
, (- s;m f ,m f ) + L^Ml) - 1 

-r f , w (Bg( - s;m f ,m f ) + L^mj,) - l) j 
- 4{ [4 + ^ [4 + (4 - r rw ) t w ] + (2 - T Jf- (4 + d^)) 
x d w R w M^Co( — m 2 ,—m 2 ,—s]M w ,mf,M w ) 
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+ \ ( 2 ~ rf ' w + [ 4 ~ r/ ' H '( 4 + d wtw)]RW^ 

x S F ( - s; M w , M w ) + L„(M£) - 1 

- H + i - 2 + t 1 + ( 4 - ^*-)^]) 



+ 2L„(M£) 



2d„ 



2 - 



while the three auxiliary functions read: 

F itm =2sC ( - mj, -m 2 f , -s; to/, A, to/) 

- 3B F ( - s; 0, 0) - 3L M (m2) + 2L A (roJ) + 2 . 



lim 



lim 



2(i + R z y 

R z 



MjC (0,0,-s;0,M z ,0) 



(3 + 2i? z ) [B f (-s;0,0)+I^(M|)-1 
2 



1 

2' 



3 + d 2 w [1 + (2 + r f . w ) R w ] 



x M w C (0,0,-s;m f ,M w ,mf) 
- \ 6 - r/' W + 2d w (2 + r/' W ) R v 



Bq(- s; m f , m f ) + L^mj,) - 1 



J- + ^+2 + 2(2 + r/^)^ v 



(2.89) 



L,(m 2 f ,) - L,(M W ) 



(2.90) 



where 



dw 



In — 

A 2 ' 

1 - rf W , 



(2.91) 



with A being the photon mass. 

The formulae, presented in this subsection, are useful to compute the one- 
loop corrections for the channels with massless or light final state fermions / 
and heavy virtual fermions /'. 
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2.6 Library of scalar form factors for electron vertex 



Besides Btt clusters, we need also Bee clusters, which can, in principle, be 
taken from j^j or derived from the formulae of previous subsection in the limit 
TO/' — ► and rrif = m e . Here we simply list the results: 



F A ' e (s) (3« e 2 + a 2 e ) a e ^*' e (*) + r- e (s) 



with 



T w ' e (s) = -^ w » fc - c (s) + 4^ H '"- e (s) . 
In Eq. I)2.92|l we use four more auxiliary functions: 



T A - e (s) 



2sC ( - TOg, to 2 , -s;TO e , A, m e ) 
35^ ( - s; 0, 0) - 3L m (to 2 ) + 2L A (m 2 ) + 2 



(s) ^ 2 (1 + i ^ } M z 2 C (0, 0, - S ; 0, M Zl 0) 
<(- S ;0,0)+L M (M 2 ) 



5"^ 



i? F (- S ;0,0)+L A1 (M 2 )-l 



r-" e (s) 



- 2i? w \Ml Co(0, 0, -s; M w , 0, M w ) 
+ -Bjf ( — s; M w , M w ) + L M (Af* ) - 1 

- C (0, 0, -s; M w , 0, M w ) 
-B£(-s; M w ,M w ) ~ 3L M (M* ) + H 



( S ) ( + m ^ Cq ( 0? o, - s; 0, M w , 0) 



Bg(-s;0,Q)+L^Ml) 



R w 



B*(-s',0,0)+L M {M*)-l 



(2.92) 
(2.93) 



(2.94) 



(2.95) 



(2.96) 



(2.97) 



where L\(m 2 e ) — ln(m 2 /A 2 ). 
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2.7 Amplitudes of boxes 

The contributions of QED AA and ZA boxes form gauge-invariant and UV- 
finite subsets. In terms of six structures (L, R) ® (L, R, D) they read: 



B L 



d-\-c 4 

= k LK ?- 

s 



[7m7+ ® 7^ 7+] Kl («,*,«) 
+ [7m7+®7*.7-]^m (*,*,«) 
+ [7^7- ® 7^7+] ( s , u ) 

+ [7/i7-®7/*7-]^M (*.*»«) ( 2 - 98 ) 
+ bV7+ ® (-iro/JD,,)] ^ (s, t, u) 
+ [7m7- ® (-imy/D^)] .Fj£ (a, t, u) 



where Llf = AA, ZA, and for shortening of presentation we factorize out 
normalization factors: 



k — s w Q e Q j , 



y,ZA _ s wQeQf 



k ZZ = 



32ct 



(2.99) 



2.7.1 AA box contribution 

There are only two A A diagrams, direct and crossed: 



AAA/WVWWWW^ 

A 




f e~ A 

Figure 10: Direct and crossed AA boxes. 



The six form factors of AA boxes might be expressed in terms of four 
auxiliary functions T\ and ^1,2,3: 

Ttt («, t, u) =T^ {a, t, u) = Hi {a, t) - Hi (a, u) + W 2 (*, t) + H 3 {a, u) , 
^ (s, t, u) =^ (s, t, u) = Hi (*, t) - Hi (s, u) - H 2 (a, u) - H 3 (a, t) , 
T$£ (a, t, u) =T A L A D (a, t, u) = T x (a, t) - T x (a, u) . (2.100) 
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The auxiliary functions are rather short: 



•MM) 



n 2 {s,t) 



n 3 (a,t) 



Here 



1 s J 1 

2 I 



t 3 _J AA (-s,-t;m e ,m f ) 



+ ts 2 C a ( - ml, -ml, -s; 0,m e , 0) 



+ t 

- 2 
+ 2 



4m^ s ( s 



2ml 



A 3r 



14,. 



Co ( — m^, —mf, —s; 0, mj, 0) 



i- 



B^(- S ;0,0)-B^(-m 2 ;m f ,0) 
B o{~ t;m e ,m f ) - Bq ( - m 2 f ;m f ,0) 



(2.101) 



Wi (M)=--r 



t + 



+ t-Co( — m l, m> 



+ 



A4 r J A4 r 

2 f ,-t;m e ,0,m f ) 



J AA (-s,-t;m e ,mf) 



(2.102) 



.1-2— 

2A 4r V i 



B(f ( - t;m e ,m f ) - Bq ( - m 2 f ;m f ,0) 



(s + 2t-) \l-t+ 



s€'o(- m 2 ., -ml, -s; 0, m e , 0) 



4A4 r I V A4 r 

2m 2 f s - (s - Am)) (^s + 2i_ - (st+ + 2tt_) -^-^j 

x Co( - m 2 , -m 2 , -s; 0,m/, 0) 
<(-s;0,0) -B^(-t;m e ,m f ) 



+ 2t 



-Am 2 



Bq( ~ t;m e ,m f ) - B^ 



m 



\;m f ,0) 



4Ai 



(s + 2t-)t 3 _J AA (—s, -t; m e , mf) . 



(2.103) 



(2.104) 



A 4r = -tu + mj, (2.105) 

and J AA (Q 2 , P 2 ; Mi, M 2 ) is due to a procedure of disentangling of the infrared 
divergences from D . Its explicit expression reads (P 2 > 0, Q 2 < 0, and M\ 
is ignored everywhere but for the arguments of In) : 
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J AA (Q 2 ,P 2 ;M 1 ,M 2 ) 



In 



In 



-Q 2 P 2 \-Q 



P 1 



2 



-Q 2 



.1/ 



+ In" ( 1 + -pf ) - 2Li 2 [ — 



P 2 



-2 hl 



" In 



M 2 



A 2 



(2.106) 



Furthermore, the relevant infrared divergent Cq function (again P 2 > 0), is 



C^ R (-M?, -M 2 2 , P 2 ; Mi, A, M 2 ) = 



1 



p2 



1 ^ 2 
2 ln 



P 2 



2Af 

im 2 





r i 


j, 


M 2 M 2 _ 



In 



A 2 



(2.107) 



P 2 



In" 1 



P 2 



-2Li 2 

where A P = P 2 + M 2 . 

2.7.2 ZA box contribution 

In Pjr gauge there are eight ZA boxes; however, since the electron mass is 
ignored, only four diagrams without 4>o contribute, see Fig. 1111 

The six relevant scalar form factors are conveniently presentable in the 
form of differences of /- and u- dependent functions: 



(2.108) 



where the index IJ is any pair of L, R®L, R, D. The 12 Fff functions depend 
on 6 auxiliary functions by means of equations where the coupling constants 
are factored out: 



•> RR \°J L 
•J~ ' RR ( S J U 

T ZA (s t 

J RD \°1 u 
3~ RD ( S J U 



-a e SfHi (s, u) + a e cr f H 2 (s, u) , 
--5 e 5 f Qi (s, t) + 5 e <7 f g 2 (s, t) , 
-6 e cTfHi (s, u) + 5 e S f 7i 2 (s, u) , 
--a e a f TLi (s,t) + a e 8/H 2 (s,t) , 
^J>fQ\ (s,u) +<T e a f g 2 (s,u), 
=S e SfHi (s, t) + 5 e a f H 2 (s, t) , 
-8 e (TfGi (s, u) + S e S f Q 2 (s,u), 
--(T e a f Ti (s,t) +<r e 5 f T 2 (s,t) , 
-<j e bfT\ (s, u) + a e a s T 2 (s, u) , 
-S^f^i (s,t) + 8 e (jfF 2 (s,t) , 
-S^fTi (s, u) + 5 e 5jT 2 (s, u) . 



(2.109) 
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Finally, we present these 6 auxiliary functions: 



Hi (s,t) 



1 s 



Fi(a,t)=-- 1 r—<t-[R z + 2 + 2t— J AZ (-s, -t; m e ,m f ) 

O A4 r y S ZA4 r J 

- (^s + + 2^-^C ( - m 2 e ,-ml-s;M z ,m e ,0) 



-2( 



tMl s-tt 
t + 2t- + - 2 



A 4r 

X 



t+ 



2^(M 2 z -4m 2 f )+s_ 1 + 2 



Co( - mj, -Wg, -i; to/, M z ,TO e ) 



^3r 



C ( - m 2 , -m 2 f , -s; M z ,m f ,0) 



2Bq( - t;m ei m f ) 

- Bg( - m); M z , to/) - B$ { - m 2 f ;m f , 0) 
2<(- S ;M z ,0) 

-B^(-m 2 f ;M z ,m f )-B^(-m 2 f ;m f ,0) 



1 s 

8A 4r 



— J A z{-s,-t;m e ,m s ) 



- 2iC ( — rrij, —m 2 , —t;mf,M z , m e 

„ m 2 



(2.110) 



7 



t - t_ i? z - 2 + 



1 



A 4r 

x J AZ (-s,—t;m e ,mf) 



- 2 



Rz 



+ 



t - t- \R Z - 2 + 



S-t-i 



1 

A 4r / 



1 



s+ + s_ (s + 2t-) 



x C (- m 2 f ,-m 2 e ,-t;m f ,M z ,m e ) 
t- 



+ 



\ 4r 

s_ - 2t + - s_ (s - 4to^) 

A 4r 



C ( - TOg, TOg , -s; M z , to 6 ,0) 



+ 2 

- 2 



x C (- m 2 ,-m 2 ,-s-,M z ,mf,0) 
Bq ( - s; M z , 0) - < ( - t; m e , m/) 

4 



2BF 



o v - t:m e ,m f ) 
BF(-m 2 ;M z ,m f )-BF(-m 2 ;m f ,0) )\, (2.111) 
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n 2 (s,t)=- 



t_ mi / 1 m 



S- 2 \ s A 4r 



J AZ (-s,-t;m e ,m f ) 



+ — C ( - ml,-mj, -t; m e ,0, m/) 



t_ TO? 



Mr 



C ( - to 2 , -m^-ijTO^M^TOe) j, 



£1 (M)= 7 



1 



s_ 2s 2A 4r 



to) + i_ 2t + - mji? z 



Mr 



x J AZ (-s,-t;m e ,mf) + — C ( - m 2 e ,-m 2 , -t;m e ,0, to/) 



— [m} + 2t-t+ + ^- 



Mr 



x C (- m 2 f ,-ml,-t;mf,M z ,rn e ) 



+ 



2A 



4r 



ts+ + - + ) C ( - m 2 , -to 2 , -s;M Zl m e , 0) 



Mr 



- (2m 2 f t. -t(a+- Am 2 ) 

x C ( - m 2 ,-m 2 ,-s;M z ,m f ,0) 
+ 2t ^Bq ( - s; M z , 0) - - i; m e , m f ) 

2.B F ( - t;m e ,m f ) 
-B${-m);M z ,m f ) - <( - m 2 ,; m/, 0)] ) }, 



£ 2 (*,*) =- 



SB!', 



1 t 

s A 4r 



J AZ (-s,-t;m e ,m f ) 



t 



- 2— Co ( - m / ; "V, M z , m e ) 

A 4r 



where a new notation were introduced for the invariants: 

^ 2 . 



s± = s ± M 2 , t±=t± m 2 f , 



and for the new function J TJ (Q 2 , P 2 ; Mi, M 2 ) ■ 



J AZ (Q 2 ,P 2 ;Mi,M 2 ) = p2 + M2 In 



Q 2 + M\ 
Ml 



In 



M 2 Mf 



(P 2 + M|) 2 



(2.112) 



(2.113) 



(2.114) 



(2.115) 
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2.7.3 The WW box 



There is only one, direct or crossed, WW diagram contributing to our process, 
see Fig. IT2l 

t 




e + w b 

9 « fVAAAAAAAAAf- 



(u, c)t 



e~ W t 

Figure 12: Crossed and direct WW boxes 



- — *wvwww\* — «- 

e- W b 



Here we give the contribution of this diagram to the scalar form factor LL: 

9 4 1 



(2.116) 



where 



(u — m b ) Dq(0, 0, — m t , — m t , — s, — u; M w , 0, M w , m&) 



+C* ( - m?, -m?, -s; M w , 0, M w ) +C* (0, 0, -s; M w , 0, M w ) 

(2.117) 

with s, £, and u being the usual Mandelstamm variables satisfying 

s + t + u = 2m 2 f , f = torb. (2.118) 
For the processes, where the direct box contributes, we have: 



>pWW _ 

11 ~~ - r (3) 



2/ 



('->- /v 2 i /_ + ( -ni; x + ) t) ^ 



^4r 



K 2 + 



Z?o(0, 0, -m^, — m&, -s, — f; M w , 0, M w , m t ) 

,4 



K 2 t+t- -Kx[^-+t 



1 + ir 1 

— + Ai— — -2- 

'4r t ^4 r 



x C (-ml,-ml,-t;M w ,m t ,M w ) 



2[K 2 t + Kit+t--^-)t- -^—Co (0, -mg, -t; 0, M w , m t ) 



43 



2 - K 2 + (K 2 U- - Kit+) t —^- + K 1 ^' ' 



t Ai 



+ 2t_ 
- 2t + 



1 

A 4r 
1 



4/e s'wQb 



t A 4r 

x C (0,0,-s;M w ,0,M w ) 
B^(-s;M w ,M w )-B^(-m 2 b ;M w ,m t ) 

B^(-t;0,m t )-B^(-m 2 b ;M w ,m t 

t'_ 1 



t 2 



(t_Y-2K 1 t+- [m z b —-K 1 
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A^ 



1 



+ Kft+t- j ^0(0, 0, -m 2 , -rog, -s, -f, M w , 0, M w , m t ) 



m t 
t 2 



i K ^ +K 4)iz- K ' 



t z + t- 1 



t Aj r 



+ 2 



X 2 + 



C ( - m 2 , -m 2 1 -s;M w ,m t ,M w ) 
— — C (0, -ra 2 , -t; 0, M w ,m t ) 

A 4r 



+ 2 



1 



1 I y M- 1 

A 4r +Al t A 2 r 



A 4r 
2t+ 1 



t_ A 



4r 



s 1 

j(3)A 4i 



x Co(0,0,-s;M w ,0,M w ) 
5<f ( - s; M w , M w ) - S F ( - to 2 ; M w , m t ) 

Bg(-t;0,m t ) - B^(-m 2 ;M w ,m t ) 



(2.120) 



2m 2 t_ + K 3 t + + K 2 t_t 



A 4r 



1 

A^ 



x D (0, 0, -m 2 , -to 2 , -s, -t; M w , 0, M w ,m t ) 
fmlt, \ 1 1 

x Co(-TOt,-TOj,-s;M w ,mt,M w ) 

r ' -1 

-* - t + -^ + 2£T 1 £_— C o (0,-mg,-t;0,M lv ,m t ) 

t t ZA47-' 

- ^3-^i^^)c (0,0,- S ;M w ,0,M w ) 



- 2t+- 



+ 2 r 



Kir 



B F (-t;0,m t ) -B^(-m 2 b ;M w ,m t ) 



(2.121) 
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where 



In the limit nib 



K 3 

t'_ 
t- 



2M' n 

9 9 

mfirif, 
t 2 

2Ml - 



t t_ 



= t 



= t — m 



6; 



4- 

f t 



= t + m t , 
= t + m 2 . 



0, for F^ L W one has: 



(2.122) 



ll lim 



where 



21, 



(3) 



57 + '|t-«i-; 



lim\2 



x D (0, 0, 0, 0, -s, -t\ M w , 0, M w , m t ) 
u it 









-)] 







Co (0, 0, -s; M w , 0, M w ) + C o (0, 0, -s; M w , m 4 , Af w ) 



1 - 



lim 



C o (0,0,-i;0,Af w ,m t ) 



m (-s;M w ,M w ) -BS(-t]0,rrH) 



(2.123) 



lim 



2M l w + i_ 



(2.124) 



2.7.4 The ZZ box 

There are four ZZ box diagrams, which form a gauge-invariant and UV finite 
subset, see Fig. IT31 

First, we use three auxiliary functions J- \Ti,Q: 



1 

J I 



s, t, u) 



RD 

-zz 

LB, 

T Z J (s,t,u) 



T (a e ,af,5f, s, t) - F(a e , 5 f ,(jf, s, u) , 

T (6 e ,6f,(Tf,S,t) - T (S e ,(Jf,Sf,S,u) , 



RL 
-ZZ 
LL 



H(a e ,af,Sf,s,t) 
H (5 e ,8f,crf, s, t) - 



- G (a e ,Sf,af,s,u) . 

G ' (5 e ,<Tf,8f,S,u) , 



£ (a e ,af,5f, s,t) — H (a e , 5f,af,s,u) . 
G (S e ,Sf, crf,s,t)-H (S e ,af, S f , s, u) . 



(2.125) 
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e + Z f e+ z f 




Figure 13: Crossed and direct ZZ boxes. 



Separating out Z fcrmion coupling constants and some common factors, 
we introduce more auxiliary functions. For T (<r e , 07, 5f, s, t) denned as 



a\a 2 T\ (s, t) + (j\o fbfT-i (s, t) 



(2.126) 



there are two 

t) = (t+Ml + t 2 _ + [(-s z + 2M 2 ) t 2 _ + UMt] i=-) 
x L> (0, 0, -m 2 f , -m 2 f , -s, -t; M z , 0, M z ,m f ) 

t- (t-t + 2M 2 z t+) + smj 2t+ (M| - 2m 2 f ) 



m}-M 2 - 



A 4r Afr 

x Cq( - m 2 ,-m 2 , s; M z , m f ,M z ) 



- 2 



2 \ t 3 + 2tt-M 2 - m 2 (3tt- + mfj 



M 2 Z - m) 



s + 2t- 



t- J ' A 4r 

x C (0,-m 2 f ,t;0,M z ,m f ) 
(t - 2t_ - 2M|) t 2 _ + sm) 



Mr 



C o (0,0,s;M z ,0,M z ) 



A 



■3r 



B Q ( - s; M z , M z ) - B ( - m 2 f ; M z ,m f ) 



2t r 



B (*;m/,0) -S (-m^;M z ,m/) 



(2.127) 



and 



^2 (a, t) = - (t 2 - + 2M|i) £> (0, 0, -m^, -m^, -s, -t; M z ,0, M z ,m f ) 
-t+C ( - m 2 f ,-m 2 f ,s;M z ,rrif,M z ) 

+ 2tC (0, -m), t; 0, M z , m f ) - t_C* o (0, 0, s; M z , 0, M z ) . (2.128) 
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For H written out as 



m 2 fS 

H(a e ,cjf,5f,s,t) =— — 

A 4r 



a 2 a 2 Hi (s, t) + a 2 afSfH2 (s, t) 
we need three auxiliary subfunctions: 



+ a 2 e 5}H 3 (s,t), 
(2.129) 



Hx (s,t) 



^--(t_ + Mlf-2Mlt- S:l1 



2A 



■4r 



x L>o(0, 0, -mj, -mj, -s, -t; M z , 0, M z ,m f ) 
. t\ + A 3r m 2 f 



t+ 



+ 2t + i_ + 



2A 4r 

2M 2 t s z t-t+ 
__ 



\ 4r 



C ( - m^, -m^, s; M z ,to/, M z ) 
C (0,0,s;M z ,0,M z ) 



t_ + 

2 



( 



s z ( sm 2 f — t 2 



2A 



4r 



Co(0,-m 2 f ,t;0,M z ,m f ) 

2^^B (t;m f ,0) - B (- m);M z ,m f )^ 

- B (t;m f ,0) + B (s;M z ,M z ), (2.130) 
H 2 (s, t) = [t- (s+t-) + 2m 2 f M 2 ] D (0, 0, -m 2 f , -m 2 f , -s, -t; M z , 0, M z ,m f ) 

-(s + t- -2m 2 f )C (-m 2 f ,-m 2 f ,s;M z ,m f ,M z ) 

— 2m 2 Ca( — mj, —m 2 , s; 0, M z ,m,f) 

+ (s + i_) C (0, 0, s; M z , 0, M z ) , (2.131) 

and 

H 3 {s,t)=-s t-D (0,0,-m 2 f ,-m 2 f ,-s,-t;M z ,0,M z ,rn f ) (2.132) 

+ C ( - m^, -m 2 f , s; M z , m f , M z ) + C (0, -m^, i; M z , 0, M z ) 
Finally, also, denned as follows: 



g(a e ,a f ,S f ,s,t) = — 



ala)Q Y (s,t)+a 2 e a f Sfg 2 (s,t) 



(2.133) 
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needs only two additional functions: 



Gi (s,t) 



t-Ul 2M 



\j2 



sit 



2A 



ir 



t- 



7 



t-t + tM 2 z {2m) - M 2 Z ) 



x £> (0, 0, -mj, -m), -s, -t; M z , 0, M z ,m f ) 
t-tl 



t+t 



Mr 



t+ {t + t-) + 2M 2 z m) — 



C ( - mj, -m 2 f ,s; M z , m f ,M z ) 
t t + t-ts z 



C (0,-m 2 f ,-t;Q,M z ,m f ) 



t-t- 



t+t 2 _ 



2mj ( 1 



A4 r 

C o (0,0,a;M z ,0,M z ) 
Bo(*;m/,Q) - B (-m 2 f ;Mz,mf) 



B (s;M z ,M z ) -B o (t;m f ,0) 



(2.134) 



and 



G 2 (s,t) = m) U 2 _ + 2M 2 ^D (Q,0,-m 2 f> -m 2 ,-s,-t;M z ,0,M z ,m f ) 

+ t+C ( - m 2 f7 -m 2 f ,s; M z ,m f ,M z ) - 2<C* (0, -m), t; 0, M z ,m f ) 



+ t-C o (0,0,s;M z ,0,M z 



(2.135) 



We recall that A4 



\i T = —tu + mff together with A 3r of Eq. I|2.84|l denotes 
remnants of Gram determinants that remained after cancellation of factors s 
and 4, leading to a simplification of the expressions. 



2.7.5 Transition to the L, Q,D basis 

Since the ZZ box contribution is given in the L, R, D basis, while all the 
rest is in the L, Q, D basis, we should transfer one of them to a chosen basis. 
At this phase of the calculations there is not much difference which basis is 
chosen. For definiteness we choose the L, Q, D basis and transfer the ZZ box 
contribution to it. The transition formulae are simple: 

TU (s, t, u) = TU is, t, u) + TU (s, t, u) - TU (*, *, «) - Kt («, *, «) > 



f z Q l (s,t,u) 



T z R l{s,t,u)-T z R z R {s,t,u) 



Tl z Q ( s > *> u ) = 2 HI ( s > *> u ) ~ Kl (s, t, u) 
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F™(s,t,u) = 2T z R z D { Sl t,u). 



(2.136) 



2.7.6 Box Born interferences 

Any box, described by the amplitude in Eq. I|2.98[l . interfering with 7 and Z 
exchange tree-level amplitudes, gives rise to two contributions to the differ- 
ential cross-sections, which are useful for internal cross-checks: 



BOX®BORN- 



cx 8Q e Q t Re 



(s + 1-) + sm 2 



{T LL + T R 

R) 



(2.137) 



+(sm} + t 2 _)(F LR + T RL )- 2m) (st + t 2 _){T LD + T RD ) 



8Re 



(s + t~) + sm 2 f 8 f {(Te^LL + 5 e T RR ) 
+ 2 (s + t-) 2 a f a e T LL + 2t 2 _a f 8 e T RL 

+ (sm 2 + t 2 _)5 f (a e T LR + 5 e J 7 RL )+2sm 2 a f (o e T LR + 5 e T RR ) 

(2.138) 



2m) (st + tt) v f (a e T LD + 5 e T RD ) X 



3 Scalar form factors for electroweak ampli- 
tude 

Having all the building blocks, it is time to construct complete electroweak 
scalar form factors. 



3.1 Vertices scalar form factors 

We begin with two vertex contributions: 





Figure 14: Electron (a) and final fermion (b) vertices in 
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In the same way as described in jS] , we reduce two vertex contributions to 
our six form factors: 



F LL (s 
F QL (s 
F LQ (a 

Fqq (S 

F LD (s 

F QD (s 
where 



Fl ee (*)+F* ff (a)-ic 2 w A{M w ), 

F z Q ee (a) + F*Jf (a) - 2c 2 w A(M w ) + k [F?" (a) - 2A(M W )] , 
Fr 00 - 2c 2 w A(M w ) + 1%" (a) + k [F2 ee (a) - 2A(M W )] , 
F z Q ee (s) + F*ff (a) - A [F2 ee (a) + F2" (a)] , 



F z J f (s), 

Fl" (a) + kF2 ff (s) , 



k = 4 {Rz - 1) ■ 

With the term containing A(M W ), 

1 M 2 
AM w ) = --ln^f 



(3.139) 



(3.140) 



(3.141) 



we explicitly show the contribution of the so-called special vertices |36| . Note 
that they accompany every L form factor. The poles 1/e originating from 
special vertices will be canceled in the sum of all contributions, including 
self- energies and boxes. 

3.2 Bosonic self-energies and bosonic counterterms 

The contributions to form factors from bosonic self-energy diagrams and 
counterterms, originating from bosonic self-energy diagrams, come from four 
classes of diagrams; their sum is depicted by a black circle in Fig. ^| Their 
contribution to the four scalar form factors is derived in [5]: 



Ft {s) =V Z (s)-s 2 w U(0) + 



Ap + Ap hos , (3.142) 



*%m (*) = v * (*) - (n^) + K°: (*) - 4^(0) - ( Ap + Ap 



-bos 



(3.143) 



F, 



ci.bos 



(a) =V h ™ (a) 2 (ll h ™(a) + n^ s (a)) + , 
- s 2 v n h ° s (0) - 4- f A P bos + Ap hos 



(o) 



(3.144) 
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Figure 15: Bosonic self-energies and counterterms for ee — > (Z, 7) — ► //. 

F%t ( s ) =^ («) ~ ™% (s) - (0) - -1 V° r - (3-145) 

We note that the term (s) — ife (s)] is conventionally extracted 

from Fqq 1ci (s). This contribution is shifted to A™ A , Eq. ((Ol) . 

In Eqs. l|3.142H3.145|> A P b ° S and n^° s (s) stand for shifts of bosonic self- 
energies. They have the same origin as special vertices and they are equal 
to: 

Ap hos = As 2 w A(M w ), (3.146) 

n^° s (s) = -2R W A{M W ), (3.147) 

see Eqs. (6.137) and (6.139) of |1 1|. These poles also cancel in the sum of 
all contributions. 
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3.3 Total scalar form factors of the one-loop amplitude 

Adding all contributions together, we observe the cancellation of all poles. 
The ultraviolet-finite results for six scalar form factors are: 

F LL (a, t, u) =TT (a) + J*" (*) + Jfz. («) + 16fe^r (*> *> «) , 
F Qi ( s , t, u) =J^- ( S ) + ^ ( s ) + fc ^ (s) + T« L (s) + 16kF£* (a, t, u) , 
F LQ (a, t, u) =Tr (a) + K SS ( s ) + k FT (*) + T? Q (a) + 16fc^°* (a, t, u) , 
F QQ (s, t, u) =^ ee (a) + (a) - A [j^ ee (s) + (s) 1 + (s) 

F i£} (s, t, u) =T"J* (a) + 16fc^° x (a, i, u) , 

F QD (a, t, u) =T Z J S (a) + k FZ ff (s) + 16kF*° x {a, t, u) . (3.148) 
For IJ = LQ, QL, QQ, LD, QD components of a box contribution we have: 
ff?* (a, t, u) =k AA T AA (a, t, u) + k ZA Ff A (a, t, u) + k*'Tff (a, t, u) 

+ k ww FZ W (s,t,u) , (3.149) 

where 

k ww = . (3.150) 

Moreover, 

E ( 3 - 151 ) 

except for J r 2 A (s) = and T2 H (s) = 0. 

4 Improved Born Approximation cross-section 

4.1 Improved Born Approximation cross-section 

In this section we give the improved Born approximation (IBA) differential 
in the scattering angle cross-section. It is derived by simple squaring the 
(■y+Z) exchange IBA amplitude, Eqs. I|1.8|l - p.l()|l . and accounting for proper 
normalization factors. We simply give the result: 

' a -P t N c (a^ A + ^ z A +o' z BA ), (4.152) 
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where (3 t = \Jl — 4mJ/s and 



ib a _ n2 I «2 

77 



(s 2 + 2at+2t_) a (a) 



(4.153) 



a^- 4 =2Q / Re<^x 2 (s + L) z + S mj F LL {s,t,u) 



+ (s 2 + 2st + 2t 2 
-4m ) (st + t 2 _) 



F QL (s, t, u) + F LQ (s, t, u) + Fqq (s, t, u) 



F LD (s,t,u)+F QD (s,t,u) 



a* (*)}, 



(4.154) 



F LL (s,t,u) +F LL (s,t,u)F*(s,t,u) 



+ 2 
+ 4 



{s + t-f+t 2 _ F QL (s,t,u) 



(s + t-) 2 + am 2 







2 




2 


F LQ {s,t,u) 


+ 


Fqq (s,t,u) 



2F LL (a,t,u)F* Q (s,t,u) 
+ F LL (s, t, u) F* Q (s, t, u) + F QL (s, t, u) F* Q (s, t, u) 
+ [s 2 + 2(st + t 2 _)} 

+ 2 (f ql (a, t, u) + F LQ (a, t, u) ) F* Q (s. i u ) 
- 8m 2 (st + t 2 _) (^2F LD (s, t, u) + F QD {a, t, u) ) F. L (».!.„) 
+ (f ld {a, t, u) + F QD {a, t, u) ^ F* L {a, t, u) 
+ (2F LD {a, t, u) + F QD (a, t, u) ^ F* LQ {a, t, u) 
+ I F LD (s, t, u) + F QD (s, t, u) ) F* Q (s, t, u) 



- 2m 2 , (at + t 2 _) A 3r 



F LD (s,t,u) 



+ 2F LD (s, t, u) F* (s, t, u) + F QD (s, t, u) 



(4.155) 
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4.2 The e + e — ► // process in the helicity amplitudes 

According to the analysis of the EW part and QED parts, we have the com- 
plete answer for the amplitude of the process e + e~ — * tt. 

The aim of this section is to adapt the helicity amplitude techniques for 
the description of our process. We produced an alternative analytic answer 
for the same amplitude using the method suggested in Ref. \A7\ . 

In general, there are 16 helicity amplitude for any 2/ — > 2/ process. For 
the unpolarized case, and when the electron mass is ignored, we are left with 
only six helicity amplitudes, which depend on kinematical variables and our 
six scalar form factors: 



A 



+++- 



=0, A+++- = 0, .4++-+ = 0, 



= 0, 



A+- 
A+- 



=s (1 - costf) ( QeQfFaa + XzS e (1 + fif) 4 3) F QL + S f F t 



QQ 



=a (1 + cos,?) ^QeQfFaa + Xz S e [(I - /3/)/ t (3) F 0£ + 5 f F QQ 
=.4 + _ ++ = 2y/smfSm'&( Q e QfF GG 



Xz$e 



4 3) F QL +8 f F QQ + -sP/irF QD 



2 r(3) ; 



A =A | = —2\fsmf sini?( Q e QjF i 



A. 



A- 



A—++ - 
Here 



+ \ z 2I^4 3) F LL + 2I^S f F LQ + S e 4 3) F QL + S e S f F c 
+ ^ sf 3 f 2 4 3) (2I^F LD +5 e F c , 
s (1 + costf) (^Q e Q f F GG + X z [ (1 + 0f) (2I^lf ] F LL + SJ^ F QL ) 
+ S f (2I^F LQ +6 e F QQ )]\ 

s (1 - cos i?) {Q e QfF GG + xz [(1 - 0f) 4 3) (2I^F LL + S e F QI }j 
+ S f (2I^F LQ +S e F QQ ^y 

0, A—+- = 0, A—+ = 0, A =0. (4.1) 



cos $ = (t — mj + 



2) si3 f 



(4.2) 
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and for the amplitude A\ i \ j \ k \ l each index \% j t k,i) takes two values (± = ±1) 
meaning 2 times the projection of spins e + ,e~,t,t onto their corresponding 
momentum. The differential cross-section for the unpolarized case is: 



da na 2 



P f N c J2 |-4a, Wi | 2 . (4.3) 



d cos $ s 3 

Xi Xj Afc A; 

We checked that this expression is analytically identical to Eq. I|4.152|) . 
Both expressions, Eq. (|4.152|l and Eq. (|4.3|) . contains, however, spurious con- 
tributions of the two-loop order (squares of one-loop terms), which should be 
suppressed, since we would like to have a complete one- loop result. 

This may be achieved with a simple trick. First of all, let us note that 
if all form factors are: F,,, = 1 for IJ = LL, LQ, QL, QQ and Fu = for 
I J = LD, QD, we have the tree level. At the one-loop level LL, LQ, QL, QQ 
form factors may be represented as: 

F„ = 1 + 5 ^F„, (4.4) 

and 

for IJ = LD, QD. 

Instead of Eq. (|4.4(l for the four form factors we write 

Fjj = Z + j^-Ftj , (4.6) 

and note that the cross section is a function of six form factors. 
Then the one-loop result is apparently equal: 

cfcr 1 ) da ,„ M . da 



l Z = 1 ]~TZZZ-ai Z = °]- ( 4 - 7 ) 



d cos i) d cos d d cos $ 

5 Annex 

5.1 QED vertices and soft-photon contributions 

Here we present virtual corrections due to QED vertices, a factorized part 
due to QED boxes and soft-photon contributions. The expressions in this 
subsection can be also cast from 

The formal structure of factorized virtual and soft contributions is as fol- 
lows: 



^virt+soft 

7T 



2 evirt-f soft . ^ ^ cvirt+soft i rvirt+soft 
2°ISR + We^f0 1F1 + W/0 FSR 

There are thus three types of contributions: ISR, FSR and IFI 



(5.1) 
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5.1.1 Initial-state radiation (ISR) 

The contributions of the initial-state QED e + e~j vertex and ISR soft photon 
contribution arc short, since the electron mass is ignored: 

mlf, J\ 1 2 3, 



= -ln-^^ie-lj --^ + -« e -2 + 4Li 2 (l) 

-i, _ ,„ I J '"r \ ,\ , 1,2 



C < = hl (— -^(l). (5.2) 



where 

/ - lr, ( __ 

i2 



^ = m ( — o ) ■ (5-3) 



5.1.2 Initial final state interference (IFI) 

This originates from contributions of QED boxes: 77, Z7 and initial-final 
state soft photons interference: 

<S£? = -21n^ln^, (5.4) 

= 21n^ln^-+ [F soh ( S ,t)-F s ° lt ( S ,u)} , (5.5) 

with 

1 o 1 o ( 2m 2 , \ 

F soft (s, t) = - -ll - - In 2 r, + 2 In r, In I 1 + j-J- \ (5.6) 

2 / n , 2 ™/\ , ,2/ *A ( st 



- ln ^ 1 + ^rJ +ln l-d +21n l-d ln l 1 + ^ 

+ 2Li 2 [ 1 - -4- V 2Li 2 ' « 



t-0+J \P+t- + 2mj 



-2Li 2 -2Li 2 (l) ; 
where we introduce the notation: 



4m 2 

0+ = 1 + 0, 
0- = 1-0, 
V = |- (5.7) 
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5.1.3 Final-state radiation (FSR) 



The contributions of one-loop QED ff-f vertex and final-state soft photon 
radiation are: 



evirt 

°FSR 



111 



A 2 

(1+^ 2 ) 



2/3 



f 4u^ 

U FSR 111 \2 



(1 + /3 2 ) 
2/3 

-I In 2 

2 

(1 + /5 2 ) 



In 77 — 1 



-/31n?7- 2 



In rj + 2 In?] In (1 - 77) + 2Li 2 (r/) + 4Li 2 (1) 
1 



(l + l 2 ) 
2/3 



2/3 
1 



In 77 — 1 







In 77 



In 2 77 + 2 In 77 In (1 - 77) + 2Li 2 (7/) - 2Li 2 (1) 



(5.8) 



The contribution of the ISR given in Eq. I|5.2|l may be obtaineded from these 
expressions in the limit mt = m e — > 0. 

5.1.4 Non-factorized final-state vertex 'anomalous' contributions 

To present this contribution let us introduce the definition 

0-1 



L n = In 



+ 1 



(5.9) 



The 'anomalous' part of the QED vertex contribution to the differential cross- 
section reads: 



da a 
dcos 7? 



4a 3 N c -lQ} [(QlQ) + 2Q e Q f v e v f Re( Xz ) + (v 2 e + a 2 ) v?\ Xz \ : 
x (st + t 2 _) Re(L„) + Q e Q f a e a f s (s + 2rj_) Re(L„Xz) 



( {vl + al) 



s (s - Am}) + 2 (st + t 2 _ 



2v e a e Vfa,fS (s + 2i_) ) |x z | 2 Re(i n ) 



(5.10) 



5.2 An alternative form of the cross-section for QED 
boxes 

Here we present some useful formulae, which are not in the main stream of our 
approach (described in the previous sections), but were used for internal cross- 
checks of the calculations of the QED part of the process under consideration. 

The QED boxes, Eqs. (|2.137|l and (|2.138|l . may be greatly simplified purely 
algebraically if the cross-sections are calculated analytically. For the sum of 
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AA and ZA boxes one may easily derive the cross-section: 



dcos-d 



2a A 

/3Q e QfN c Rc i QlQ'fJ 7 , 



QeQfXz 



(5.11) 



v e v f (T* + Ti v ) + a e a f (T* A + Q v ) 



\Xz 



{v\ + a 2 e ) ( vjH v + a 2 fH A ) + 2a e v e a f Vf (Q v + Q A ) 



where Xzi s ) is defined by Eq. i|1.12[l and the six cross-section form factors 
are: 



with 



T v it) = 



Tv = 


=T V (t) 


- T v (u) , 


Ta = 


--Ta (t) 


f T A (u) > 


ri v = 


--n v (t) 


- n v (u) , 


TL A = 


--Ha (t) 


- H A (u) , 


Q v -- 


--Q v (t) - 


f Qv (u) , 


Qa -- 


--Qa (t) - 


r Qa («) , 


+ (s 


+ 2t_) 


J AA {-S, - 


m 2 e ,- 


m 2 el —s 


0,m e ,0) 



-<— 2m' 
s\ 4 

si 

+ 2 ( 

's-Aml 4m?\, 



2tm^ 
2 r 



2 ' A 3r 

F I 2 



B^(-m};mf,0)-B^(-s;0,0) 



sm 



t- 



Bq( - t;m e ,mf) - -Bjf(- m^;m/,0) 



(* + *-) 



^ (t) =- 



1 f s + 2t_ 



B^(~t;m e ,m f ) -<(- S ;0,0) 
t-J AA {-s,-t;m e ,m,f) 
Q sC °( - -m^ , -s; 0, m/, 0) 

.Bjf (- t;m e ,mf) - -Bjf( - m 2 f ; m/, 0) 
Bg(-t;m e ,m f ) - - a; 0, 0) 



r + 1 

(«+*-) 



(5.12) 



(5.13) 



(5.14) 
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Ho (t) (t- + ( s + £-) 2 ) ^-JAA{-s,-t;m e ,mf) 



+ t. 



J AZ (-s,-t;m e ,mf) 



+ 



Mi + 1. 



s_ (s+ + 2t)- (ti + (s + t_) 2 )]^ 

x C (- m^, -m 2 e ,-t;m f ,M z ,m e ) 

+ s_ ^si C ( - m 2 , -ml, -s; M z ,m e ,0) 

+ stC ( - m 2 f ,-m 2 f , -s;M z ,m f ,0) 

2B o( - t;m e ,m f ) - Bq(- m 2 ;m f ,0) 

-B^(-m 2 ;M z ,m f ) 



sm f 



-(« + *-) B^{-t;m e ,m } ) - Bg (- s;M z ,0) 



(5.15) 



H v (t) =W (t) - 



2t_m 2 



J AZ (-s,-t;m e ,mf) - -J AA (-s,-t;m e ,m f ) 



+ C (- m 2 ,-ml,-t;m f ,M z ,m e ) 



2m 2 s-t 



1 - ^-JCo( - m 2 f ,-m 2 f ,-s;M z ,rrif,0) 



*3r 



2B F (-s;M z ,0) -B^-m^m/.O) 



Bl{-m 2 ;M z ,m f ) 



(5.16) 



2ml , 
W„ (t) -Ho + 



M 2 J AZ (-s,-t;m e ,m f ) 



s-J AA (-s,-t;m e ,m f ) 



+ (ss+ + t-s)C ( - m 2 , -ml, -t;m f , M z , m e ) 



-tC ( 



m f ,-m f 



,-s;M z ,mf,6) 



Bg(-t;m e ,m f ) + B* (- s; M z ,0) 



(5.17) 
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Qv (t) = 1 - ^ (s + 2t_) t-J AA (-s, —t; m e , m/) 



s + 2t_ - — (s+ + 2t_) 
is 



J AZ (-s,-t;m e ,m f ) 



S-t. 

s 



(s + 2i_) t_ - (s + + 2t_) - s_TO/ ( 2 + 



M 2 



S 1 2 M ' 



x C (-my,-me, -t;m/,Mz,m e ) 
C ( - me,— me,— s; M z ,m e ,0) 

+ C ( - trip -mj, -s; M z ,to/,0) 



- 2B ( - t;m e ,m f ) + Bq ( - m 2 ; m/, 0) 
+ Bf (-m^M^m/) 
' " ' ' [b o f ( - t; me, m/) - B F ( - s; M*, 0) 



(5.18) 



G A (t)=g v (t)-m 2 f — 2s-C (-m 2 f ,-m 2 e ,-t-,rn f ,M z ,m e ) (5.19) 
+ s+C ( - m 2 , -m 2 , -s; M z , to/, 0) 

+ 2B^(~t;m ei m f ) - ( - m 2 f ;M z ,m f ) - B$( - m 2 f ; to/, 0) 

Note the cancellation of A4,., leading to a great simplification of the expres- 
sions. Equations H5.11ll5~T^| were coded as a separate branch of eeffLib 
and, together with the vertex QED contributions described in the previous 
subsections, were used for an internal cross-check of the QED part of the 
calculations. 

Some factorized part of the AA and ZA boxes contribution is not included 
in Eq. i|5.11|) . It has the form 



-—QeQfd 



d cos ■& 7T 



virt 

IFI ' 



(5.20) 



where 5™* is given by Eq. (|5.4|l . 

The whole QED contribution can be written as follows 



da QE 



da BOT 



rvirt+soft 



+ 



da a da E 



+ 



dcosfl dcosi? dcosi) dcosi) 

where ^irt+soft ig d e fi nec j by Eqs. f5~T|) - (|5^|) . 



(5.21) 
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6 Numerical results and discussion 



All the formulae derived in this article are realized in the FORTRAN code with 
the tentative name eef f Lib. Numbers presented in this section are produced 
with the February 2002 version of the code. As compared to the December 
2000 version, used to produce numbers for Ref. |23| . the current version con- 
tains full QED corrections together with the soft-photon contribution to the 
angular distribution da / dcosfi. In this section we present several examples 
of numerical results. In particular, we will show a comparison of the elec- 
troweak form factors (EWFF) including QED corrections between eeffLib 
and another FORTRAN code, which was automatically generated from form log 
files with the aid of system s2n.f (symbols to numbers), producing a FORTRAN 
source code — a part of our SANC system. This comparison provides a powerful 
internal cross-check of our numerics that practically excludes the appearance 
of bugs in numerical results. 

We begin with showing several examples of comparison with ZFITTER 
v6.30 (2H1- In the present realization, eeffLib does not calculate M w from 
fj, decay and does not precompute either Sirlin's parameter Ar or the total 
Z width, which enters the Z boson propagator. For this reason, the three 
parameters: M w , Ar , T z were being taken from ZFITTER and used as INPUT 
for eeffLib. Moreover, present eeffLib is a purely one-loop code, while it 
was not foreseen in ZFITTER, to access just one-loop form factors with the aid 
of users flags. To accomplish the goals of comparison at the one-loop level, 
we had to modify the DIZET electroweak library. The most important change 
was an addition to the SUBROUTINE R0KANC: 

* 

* For eett 
* 

FLL=(XR0K(1)-1D0+DR )*R1/AL4PI 
FQL=FLL+(XR0K(2)-1D0)*R1/AL4PI 
FLQ=FLL+(XR0K(3)-1D0)*R1/AL4PI 
FQQ=FLL+ (XR0K (4) -1D0) *R1/AL4PI 

with the aid of which we reconstruct four scalar form factors from ZFITTER's 
effective couplings p and k's (F ld and F QD do not contribute in the massless 
approximation) . 

6.1 Flags of eeffLib 

Here we give a description of flags (user options) of eeffLib. While creating 
the code, we followed the principle to preserve as often as possible the meaning 
of flags as described in the ZFITTER description 5,. In the list below, the 
comment 'as in ZFD' means that the flag has exactly the same meaning as 
in Here we describe the set of flags of the February 2002 version of 

eeffLib. 



Gl 



ALEM=3 ! as in ZFD 
ALE2=3 ! as in ZFD 

VPOL=0 ! =0 a(0); =1,=2 as in ZFD; =3 is reserved for later use 
Note that the flag is extended to VP0L=0 to allow calculations 'without 
running of a'. 
QCDC=0 ! as in ZFD 

IT0P=1 ! as in DIZET (internal flag) 
GAMS=1 ! as in ZFD 

WEAK=1 ! as in ZFD (use WEAK=2 in v6.30 to drop out some higher- 
order terms) 

IM0MS=1 ! =0 a-scheme; —1 GFermi-scheme 

New meaning of an old flag: switches between two renormalization 

schemes. 

B0XD=6 

Together with WEAK=0 is used for an internal comparison of separate 
boxes and QED contributions: 
BOXD ! =1 with A A boxes 

! =2 with ZA boxes 

! =3 with AA and ZA boxes 

! =4 with all QED contributions 
Together with WEAK=1 (working option), it has somewhat different mean- 
ings: 

BOXD ! =0 without any boxes 
! =1 with A A boxes 
! =2 with ZA boxes 
! =3 with AA and ZA boxes 
! =4 with WW boxes 
! =5 with WW and ZZ boxes 
! =6 with all QED and EW boxes 

'Treatment' options. 
GAMZTR=1 treatment of T z . 

The option is implemented for the sake of comparison with FeynArts: 
GAMZTR=0 T z = 
GAMZTR=1 T z ^ 

EWFFTR=0 treatment of EW form factors. 

Switches between form factors and effective ZFITTER couplings p and 
k's. The option is implemented for comparison with ZFITTER: 
EWFFTR=0 electroweak form factors 
EWFFTR=1 effective couplings p and k's 
FERMTR=1 treatment of fermionic masses 

Switches between three different sets of 'effective quark masses': 
FERMTR=1 a 'standard' set of fermions masses 
FERMTR=2 , 3 'modified' sets 
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VP0LTR=1 treatment of photonic vacuum polarization 

Switches between lowest-order expression a(s) = a [1 + Aa(s)] and its 

'resummed' version, a(s) = a/ [1 — Aa(s)]: 

VP0LTR=0 lowest order 

VP0LTR=1 resummed 

EWRCTR=2 treatment of electroweak radiative corrections 
Switches between three variants for vertex corrections: 
EWRCTR=0 electroweak form factors contain only QED additions 
EWRCTR=1 electroweak form factors do not contain QED additions 
EWRCTR=2 electroweak form factors contain both QED and EW addi- 
tions 

EMASTR=0 treatment of terms with ln(s/rag) in AA and ZA boxes, which 
are present in various functions but cancel in sum: 
EMASTR=0 these terms are suppressed in all functions they enter 
EMASTR=1 these terms are retained in all functions, which results in 
loss of computer precision owing to numerical cancellation; results for 
EMASTR=0 and EMASTR=1 are equal 

EWWFFV=1 treatment of vertex and box diagrams with virtual W boson, 
switches between two variants: 

EWWFFV=0 variant of formulae without 6-quark mass 
EWWFFV=1 variant of formulae with finite 6-quark mass 

Options affecting QED contributions. 

variants of inclusion virtual and soft photon QED contribu- 

only initial-state radiation (ISR) 
only initial-final interference (IFI) 
only final-state radiation (FSR) 
all QED contributions are included 

is active only if IQED=2 or 4 and affects only Eq. (5.11): 
AA boxes interfering with 7 exchange BORN 
AA boxes 
ZA boxes 

or 4 AA+ZA boxes 



IQED= 


A 


tions: 




IQED= 


=1 


IQED= 


■2 


IQED= 


■3 


IQED= 


A 


IB0X= 


A 


IB0X= 


=0 


IB0X= 


=1 


IB0X= 


■2 


IB0X= 


■3 



6.2 eef f Lib— ZFITTER comparison of scalar form factors 

First of all we discuss the results of a computation of the complete EW part 
of the four scalar form factors (i.e. with WW and ZZ boxes), 

F LL (s,t), F QL (s,t), F LQ (s,t), F QQ (s,t), (6.1) 

for three channels: e+e~ — ► uu, dd and bb for light final fermion masses (we 
set m u — m,d = 0.1 GeV) and for bb channel we use the formulae in the limit 
mf — 0). We remind that ZFITTER is able to deliver only massless results. 
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In this comparison we use flags as in subsection 16. II and, moreover, 

M w = 80.4514958 GeV, 
Ar = 0.0284190602 , 

T z = 2.499776 GeV. (6.2) 

The form factors are shown as complex numbers for the three c.m.s. energies 
(for t = mf — s/2) and for /i = M w 2 . Tabled shows very good agreement 
with ZFITTER results (up to 6 or 7 digits agree). It should be stressed that 
total agreement with ZFITTER is not expected because in the eeffLib for 
e + e~ — > uu and dd channels we use massive expressions to compute the 
nearly massless case. Certain numerical cancellations leading to losing some 
numerical precision are expected. We should conclude that the agreement is 
very good and demonstrates that our formulae have the correct m% — * limit. 
Note, that there is an all digits agreement for e + e~ — ► bb channel since in 
both cases one uses the formulae in the limit nrj = 0. 

6.3 eeffLib— ZFITTER comparison of IBA cross-section 

As the next step of the comparison of eeffLib with calculations from the 
literature, we present a comparison of the IBA cross-section. 

In Table [3 we show the differential cross-section Eq. ()4. 15211 in pb for 
four leptonic channels and in Table El for three quarkonic channes (bb channel 
being shown twice for massless and massive b quarks) for three values of 
cos-i? = —0.9, 0, +0.9 with running e.m. coupling a (s). Since the flag setting 
VP0L=1, which is relevant to this case, affects the ZFITTER numbers, we now 
use, instead of Eq. I|6.2fl . the new INPUT set: 

M w = 80.4467671 GeV, 
Ar = 0.0284495385 , 

T z = 2.499538 GeV. (6.3) 

The numbers shown in first two rows of Tabs. H and |3] exhibit a very good 
level of agreement with ZFITTER. 

Finally, in Table 0] we give a comparison of the cross-section integrated 
within the angular interval | cos i9| < 0.999. (Flag setting is the same as for 
Table 4.) 

A typical deviation between eeffLib and ZFITTER is of the order ~ 10~ 6 , 
i.e. of the order of the required precision of the numerical integration over 
cos Examples of numbers obtained with eeffLib, which were shown in this 
section, demonstrate that ZFITTER numbers are recovered for light m t . 



2 In the preprints 1231 and 1241 we presented this and the following tables for the three 
values of scale fi = M w /10, M w , 10M W and demonstrated the scale independence of 
eeffLib numbers. 
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Table 1: EWFF for the process e+e — » //. eef f Lib-ZFITTER comparison. 





100 GeV 


200 GeV 


300 GeV 


FF 








uu channel, m u = 0.1 GeV, ma = 


ZF 


12.89583 - il.84786 
12.89583 - zl.84786 


8.24737 -ilO.64653 
8.24736 -ilO.64651 


8.98371 -il2.88466 
8.98370 -il2.88466 


Fql 

ZF 


29.30446 + i3.67330 
29.30445 + i3.67330 


29.38217 + i2. 27613 
29.38216 + i2. 27613 


31.59711 + il. 59304 
31.59710 + il. 59304 


Flq 
ZF 


29.10831 + i3. 26973 

29.10832 + i3. 26973 


29.48511 + i0.92311 

29.48512 + i0.92312 


31.65835 - i0.89711 
31.65835 - i0.89711 


Fqq 
ZF 


44.88228 + i8.85688 
44.88228 + i8.85688 


43.31854 + i9.48286 
43.31854 + i9.48286 


44.18773 +il0. 25197 
44.18773 +il0. 25196 


drf channel, = 0.1 GeV, m„ = 


F LL 
ZF 


13.70781 - il. 51002 
13.70781 - il. 51002 


15.18630 - i3. 93706 
15.18629 -i 3.93706 


8.86000 - il.80409 
8.86000 - il.80409 


Fql 
ZF 


29.64340 + i4.12394 
29.64340 + i4.12394 


31.96819 + i6. 97877 
31.96818 + i6. 97877 


31.69945 + i8.03876 
31.69944 + i8.03876 


Flq 
ZF 


29.12112 + i3. 22780 

29.12113 + i3. 22780 


30.28990 + il. 73736 
30.28990 + il. 73733 


31.69321 - i0.07000 
31.69323 - i0.07001 


Fqq 
ZF 


44.87608 + i8. 79014 

44.87609 + i8. 79014 


43.94906 + i9. 86523 
43.94905 + i9.48286 


44.19463 + ilO.38253 
44.19462 +il0.38251 


66 channel, — 0, m t — 173.8 GeV 


ZF 


11.16367 - i0.65244 
11.16367 - i0.65244 


14.68726 - il. 85834 
14.68726 - il. 85834 


11.26568 - i3.38026 
11.26568 - i3.38026 


Fql 
ZF 


26.88634 + i4.76865 
26.88634 + i4.76865 


28.08780 + i5. 27877 
28.08780 + i5. 27877 


31.02060 + i3.91962 
31.02060 + i3.91962 


Flq 
ZF 


29.12113 + i3. 22780 
29.12113 + i3. 22780 


30.28990 + il. 73733 
30.28990 + il. 73733 


31.69323 - i0.07001 
31.69323 - i0.07001 


ZF 


44.87609 + i8. 79014 
44.87609 + i8. 79014 


43.94905 + i9. 86523 
43.94905 + i9. 86523 


44.19462 + ilO.38251 
44.19462 +il0.38251 
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Table 2: Comparison of the differential EW cross-sections, [pb], leptonic chan- 
nels. First row - ZFITTER, second row - eef f Lib. 





100 GeV 


200 GeV 


300 GeV 


400 GeV 


700 GeV 


1000 GeV 


cos?? 














vv channel, m v — 


-0.9 


49.100086 
49.100085 


0.579630 
0.579630 


0.198147 
0.198147 


0.101625 
0.101625 


0.030297 
0.030297 


0.014381 
0.014381 





31.834491 
31.834490 


0.358941 
0.358941 


0.122362 
0.122362 


0.062455 
0.062455 


0.018431 
0.018431 


0.008648 
0.008648 


0.9 


66.267213 
66.267213 


0.739640 
0.739640 


0.234000 
0.234000 


0.115572 
0.115572 


0.032285 
0.032285 


0.014601 
0.014601 


e e + channel, m e = 


-0.9 


7.991109 
7.991110 


0.533194 
0.533194 


0.286463 
0.286463 


0.169610 
0.169609 


0.058032 
0.058032 


0.028672 
0.028672 





19.773604 
19.773604 


1.121915 
1.121915 


0.474004 
0.474004 


0.262276 
0.262276 


0.084408 
0.084408 


0.040933 
0.040933 


0.9 


63.501753 
63.501750 


3.504349 
3.504348 


1.466665 
1.466665 


0.811430 
0.811430 


0.263378 
0.263378 


0.128882 
0.128882 


channel, = 0.106 GeV 


-0.9 


7.991040 
7.991028 


0.533195 
0.533195 


0.286463 
0.286463 


0.169610 
0.169610 


0.058032 
0.058032 


0.028672 
0.028672 





19.773498 
19.773506 


1.121915 
1.121915 


0.474004 
0.474004 


0.262276 
0.262276 


0.084408 
0.084408 


0.040933 
0.040933 


0.9 


63.501436 
63.501422 


3.504346 
3.504345 


1.466665 
1.466664 


0.811430 
0.811430 


0.263378 
0.263378 


0.128882 
0.128882 


t + t channel, m T = 1.77705 GeV 




100 GeV 


200 GeV 


300 GeV 


400 GeV 


700 GeV 


1000 GeV 


-0.9 


7.971611 
7.968192 


0.533509 
0.533295 


0.286519 
0.286477 


0.169627 
0.169613 


0.058034 
0.058032 


0.028673 
0.028672 





19.743430 
19.745968 


1.121827 
1.121978 


0.473992 
0.474021 


0.262273 
0.262282 


0.084407 
0.084408 


0.040933 
0.040933 


0.9 


63.412131 
63.408973 


3.503720 
3.503524 


1.466558 
1.466520 


0.811398 
0.811385 


0.263375 
0.263373 


0.128882 
0.128881 
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Table 3: Comparison of the differential EW cross-sections, [pb], quarkonic 
channels. First row - ZFITTER, second row - eef f Lib. 





fOO GeV 


200 GeV 


300 GeV 


400 GeV 


700 GeV 


1000 GeV 


cos?? 














uu channel, m u = O.f GeV 


-0.9 


45.404742 
45.404602 


0.386966 
0.386966 


0.225923 
0.225923 


0.138065 
0.138065 


0.048621 
0.048621 


0.024156 
0.024156 





60.382423 
60.382566 


1.882835 
1.882837 


0.771939 
0.771939 


0.421410 
0.421410 


0.133475 
0.133475 


0.064245 
0.064245 


0.9 


f73.4675f7 
f73.46755f 


6.450000 
6.450000 


2.510881 
2.510881 


1.346620 
1.346620 


0.417295 
0.417295 


0.198842 
0.198842 


dd channel, md = 0.1 GeV 


-0.9 


86.5544f4 
86.554ff0 


0.430807 
0.430807 


0.136720 
0.136720 


0.069644 
0.069644 


0.020899 
0.020899 


0.009978 
0.009978 





72.820806 
72.820908 


1.180211 
1.180212 


0.419410 
0.419411 


0.219070 
0.219070 


0.066312 
0.066312 


0.031351 
0.031351 


0.9 


f76.7f7336 
f76.7f7376 


3.770861 
3.770861 


1.469194 
1.469194 


0.796291 
0.796291 


0.255118 
0.255118 


0.124714 
0.124714 


dd channel, mt = O.f GeV 


-0.9 


84.97033f 
84.970030 


0.416155 
0.416155 


0.136467 
0.136467 


0.071144 
0.071144 


0.020786 
0.020786 


0.009868 
0.009868 





71. 762397 
71.762496 


1.185284 
1.185285 


0.426009 
0.426009 


0.219777 
0.219777 


0.067063 
0.067063 


0.031698 
0.031698 


0.9 


174.716111 
174.716145 


3.819650 
3.819651 


1.442852 
1.442852 


0.764945 
0.764945 


0.249501 
0.249501 


0.123086 
0.123086 


dd channel, nib = 4.7 GeV 


-0.9 


84.350009 
84.084519 


0.416665 
0.415686 


0.136590 
0.136435 


0.071186 
0.071137 


0.020792 
0.020786 


0.009870 
0.009868 





71.232935 
71.248102 


1.183490 
1.184186 


0.425748 
0.425882 


0.219707 
0.219745 


0.067058 
0.067061 


0.031697 
0.031697 


0.9 


173.420486 
173.041374 


3.812786 
3.811860 


1.441766 
1.441592 


0.764642 
0.764577 


0.249478 
0.249464 


0.123084 
0.123077 
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Table 4: Comparison of the total EW cross-sections, [pb]. First row - 
ZFITTER, second row - eef f Lib. 



100 GeV 


200 GeV 


300 GeV 


Otot 


Cfb 


Ctot 


Cfb 


Ctot 




vv channel, rriy = 


84.81710 


9.509864 


0.963362 


0.089665 


0.320985 


0.021592 


84.81710 


9.509865 


0.963362 


0.089665 


0.320985 


0.021592 


e e + channel, m e = 


52.61662 


30.78899 


2.980668 


1.654673 


1.276008 


0.648414 


52.61662 


30.78899 


2.980667 


1.654673 


1.276008 


0.648414 






channel, m M = 0.106 


GeV 


52.61634 


30.78885 


2.980667 


1.654671 


1.276008 


0.648414 


52.61634 


30.78885 


2.980667 


1.654671 


1.276008 


0.648414 




T + r - 


channel, m 


T = 1.77705 GeV 




52.53632 


30.75010 


2.980435 


1.654149 


1.275972 


0.648324 


52.53661 


30.75024 


2.980443 


1.654156 


1.275974 


0.648326 


uu channel, m u = 0.1 GeV 


160.8980 


70.98406 


5.021808 


3.360848 


2.031754 


1.269556 


160.8981 


70.98416 


5.021810 


3.360848 


2.031754 


1.269556 


dd channel, m<j = 0.1 GeV 


193.7658 


50.03208 


3.120724 


1.867871 


1.149479 


0.725581 


193.7658 


50.03227 


3.120725 


1.867871 


1.149479 


0.725581 


bb channel, mj, = 0.1 GeV 


191.0416 


49.76543 


3.134547 


1.892628 


1.149243 


0.720476 


191.0416 


49.76562 


3.134547 


1.892629 


1.149243 


0.720476 


bb channel, = 4.7 GeV 


189.6321 


49.39098 


3.129824 


1.888530 


1.148541 


0.719805 


189.3855 


49.32791 


3.129858 


1.888542 


1.148565 


0.719802 



We conclude this subsection with a comment about the technical precision 
of our calculations. We do not use the looptools package [HHI- For all PV 
functions but one, namely the Dq function, we use our own coding, where 
we can control precision internally and, typically, we can guarantee 11 digits 
precision. For the Dq function we use, instead, REAL*16 TDPAZ0 coding [S] 
and the only way that is accessible to us to control the precision is to compare 
our results with those computed with the looptools package. This was done 
for typical Do functions entering the ZZ box contributions. We obtained an 
agreement within 14 or 15 digits between these two versions for all y/s = 400- 
10000 GeV and costf = 0.99, 0, -0.99. 
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6.4 Comparison with a code generated by s2n.f 

Here we present a numerical comparison of the complete scalar form factors of 
Eq. (|3.148|) . extracted from two independently created codes: the 'manually 
written' eeffLib and a code 'automatically generated' by the s2n.f software. 
We use a special input parameter set here: all lepton masses a and a conver- 
sion factor from GeV -1 to pb are taken from the 2000 Particle Data Tables 
while for quark and photon and gauge boson masses we use: 

m u , d>c , s ,j,6 = 0.062, 0.083, 1.50, 0.215, 173.8, 4.70 GeV, 

A = 1 GeV, M z = 91.1867 GeV, M w = 80.4514958 GeV. (6.4) 

As seen from Table El the numbers agree within 11-13 digits, i.e. REAL*8 
computational precision is saturated. The form factors F LD QD are multiplied 
by 10 4 to make more digits visible. 



Table 5: EWFF for the process e + e 
- s2n.f . 



tt. First row - eeffLib, second row 



7j 

j i9 | F~ 



F QL 
F LQ 

F QQ 

F LO 
F QO 



68.36399900074 - 
68.36399900068 - 
75.12465846647 + 
75.12465846641 4- 
81.01546270426 + 
81.01546270420 4- 
225.63977621858 + 
225.63977621832 + 
-0.57522852857 + 
-0.57522852857 4- 
0.16677424366 - 
0.16677424366 - 



il. 24743850729 
il. 24743850728 
34.81991916400 
34.81991916400 
19.81343626967 
19.81343626968 
154.37838168488 
154.37838168491 
iO. 34010611241 
iO. 34010611241 
0.34326069364 
0.34326069364 



79.63957322115 - 
79.63957322113 - 
76.19283172015 + 
76.19283172013 + 
82.67283873006 + 
82.67283873004 4- 
207.09189805263 + 
207.09189805254 4- 
-0.33030593699 + 
-0.33030593699 + 
0.29925308488 - 
0.29925308488 - 



20.53758995637 

20.53758995637 

28.44336684106 

28.44336684106 

13.79952080171 

13.79952080171 

133.45188150116 

133.45188150117 

0.14897150833 

0.14897150833 

0.14107543098 

0.14107543098 



F LL 
F QL 
F LQ 

F QQ 

F L D 
F QD 



48.42950001713 + 
48.42950001707 + 
68.02678564355 + '. 
68.02678564349 + '. 
73.37133716227 4- i 
73.37133716220 + i 
196.60425612149 4- -, 
196.60425612123 4- 
-0.56319765502 + 
-0.56319765502 + 
0.15893936555 - I 
0.15893936555 - I 



iS. 26103890366 

iB. 26103890367 

37.08805801477 

37.08805801477 

22.69397728402 

22.69397728403 

162.74818773960 

162.74818773963 

0.33645326768 

0.33645326768 

0.37254018572 

0.37254018572 



28.23570422021 4- 
28.23570422019 + 
58.00469565609 + i 
58.00469565607 + ', 
62.40775508619 + i 
62.40775508616 + i 
132.63279537966 + '. 
132.63279537957 + I 
-0.29067043403 + 
-0.29067043403 + 
0.26429138671 - ■ 
0.26429138671 - 



i2.43705570966 

il. 43705570966 

33.82433896562 

33.82433896562 

20.75544388763 

20.75544388764 

152.68259938740 

152.68259938741 

0.13992893252 

0.13992893252 

0.15437851127 

0.15437851127 



F QL 
F LQ 

F QQ 

F L D 
F QD 



35.17736865724 + 
35.17736865718 + 
61.03099608330 + 
61.03099608324 + 
66.08215572935 + 
66.08215572929 + 
167.63393504156 + 
167.63393504130 + 
-0.56772633347 + 
-0.56772633347 + 
0.18031346246 - 
0.18031346246 - 



14.84038724783 

14.84038724784 

39.09196533610 

39.09196533611 

25.04151178684 

25.04151178685 

170.36384103672 

170.36384103675 

:0. 34299744419 

:0. 34299744419 

0.40091423652 

0.40091423652 



0.21531292996 + 
0.21531292994 + 
40.77942026097 + 
40.77942026095 + 
44.50915974057 + 
44.50915974055 + 
59.87568281297 + 
59.87568281288 + 
-0.32035310873 + 
-0.32035310873 + 
0.34968026058 - 
0.34968026058 - 



13.66645015866 

13.66645015866 

37.94118444135 

37.94118444136 

25.51875704261 

25.51875704261 

168.13599380718 

168.13599380719 

iO. 14419510235 

■;0. 14419510235 

.0.16945266925 

.0.16945266925 
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The next Table shows a comparison of eef f Lib-s2n.f for the com- 
plete one-loop differential cross-sections da^ 1 / d cos ?? , for the standard input 
parameter set Eq. (|(j.4fl . As seen, numbers agree within 12 or 13 digits. 



Table 6: ^ os ^ f° r the process e + e — > tt. eef fLib-s2n.f comparison. 





400 GeV 


700 GeV 


1000 GeV 


COS?? 








-0.9 


0.22357662754774 
0.22357662754769 


0.06610825350063 
0.06610825350063 


0.02926006442715 
0.02926006442715 


0.0 


0.34494634728716 
0.34494634728707 


0.14342802645636 
0.14342802645634 


0.06752160108814 
0.06752160108813 


0.9 


0.54806778978208 
0.54806778978194 


0.33837133344667 
0.33837133344664 


0.16973989931024 
0.16973989931023 



6.5 About a comparison with the other codes 

As is well known, the one-loop differential cross-section of e + e~ — ► tt may be 
generated with the aid of the FeynArts system (22] ■ Previous attempts to com- 
pare with FeynArts are described in |23j. In December 2001, we were provided 
with the numbers computed with the FeynArts system |41j for da/d cos??, 
with and without QED contributions, at ^fs = 700 GeV and three values 
of cos?? = 0.9, 0, —0.9. With the current versions of eef f Lib and s2n.f 
we have 11 digits agreement for both the tree-level and one-loop-corrected 
cross-sections. 



Table 7: f a 9 for the process e+e~ -> tt with soft photons, £™ ax = Vs/10. 





400 GeV 


700 GeV 


1000 GeV 


COS?? 








-0.9 


0.17613018248935 


0.05199100267864 


0.02310170508071 


-0.5 


0.21014509428358 


0.06560630503586 


0.02882301902010 


0.0 


0.27268108572063 


0.11496514450150 


0.05495088904853 


0.5 


0.35592722356682 


0.19615154401629 


0.09941700898317 


0.9 


0.43637377538440 


0.27915043976042 


0.14426233253975 



Recently, the Bielefeld-Zeuthen team [331 performed an alternative cal- 
culation using the DIANA system Working in close contact with this 
team, we managed to perform several high-precision comparisons, reaching 
for separate contributions an agreement up to 10 digits. 

The results of a comparison between FeynArts and the Bielefeld Zcuthen 
team are presented in detail in |34|. 

As another example we present in Table0the same cross-section [ d ^ s $ ] 
as given in the tables of [32]. F° r the complete cross-section, including soft 
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photons, we agree with the Bielefeld Zcuthen calculations within 8 digits. 
(See also j3SJ.) 

Two more graphical examples of the differential and the total EWRC for 
the process e + e~ — > ti are shown in the figures and IT71 correspondingly. 
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ArWfs fMdt 

Figure 16: Relative EWRC 5{sTs, cos-d) = , /n A,, 1 to the e+e" -» ti 

da^ } >{s)/dt 

differential cross-section. Numbers near the curves show y/s in GeV. 
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<5[%] 




Figure 17: Relative EWRC to e+e" -» for M„ = 100 GeV (solid line) and 
M H = 1000 GeV (dashed line). 
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